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It  NOTES 

1  ,  SR  ON  SO  WING  Mi  t  ’  A  *V  *  A(  ti  VI  T  • 

Naval  Ordnance  Systems  C  ommand 

Washington.  I).  (  .  20360 

In  situ  measurements  of  congressional  (sound!  velocity  and  attenuation  were 
made  in  the  sea  floor  off  San  Diego  in  water  depths  between  4  and  1  100  m:  fre¬ 
quencies  were  between  3.5  and  100  kll/,  Sediment  types  ranged  from  coarse  sand 
to  elayey  silt.  These  measurements,  and  others  from  the  literature,  allowed  analyses 
of  the  relationship  between  attenuation  and  frequency,  relationships  between  attenua¬ 
tion  and  other  physical  properties,  and  study  of  appropriate  viscoelastic  models  which 
can  be  applied  to  saturated  sediments.  Some  conclusions  are:  (1)  attenuation  in 
dB/unit  length  is  approximately  dependent  on  the  first  power  of  frequency:  (2!  velocity 
dispersion  is  negligible,  or  absent,  in  water-saturated  sediments:  (3)  intergrain 
friction  appears  to  be.  by  far.  the  dominant  cause  of  wave-energy  damping  in  marine 
sediments;  viscous  losses  due  to  relative  movement  of  pore  water  and  mineral  structure 
are  probably  negligible:  (4)  a  particular  viscoelastic  model  (and  concomitant  equations) 
is  recommended:  the  model  appears  to  apply  to  both  water-saturated  rocks  and  sedi¬ 
ments:  and  (5)  a  method  is  derived  which  allows  prediction  of  compressioual-wavc 
attenuation,  given  sediment  mean  grain  si/e  or  porosity. 
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SUMMARY 


PKOULI-At 

Determine  .mil  ■%( tul \  those  characteristics  or  the  sea  floor  that  affect  the  piopaii.i 
lion  of  acoustic  energy  in  the  sea.  Support  underwater  acoustics  e\pet intents  and  theory 
by  furnishing  information  on  the  mass  physical  properties  ol  sea-floot  sediments  spe 
cil'ically.  study  the  relationships  between  tivapiency  and  sound  atleuualion.  furnish  values 
of  attenuation  in  specific  sediment  types.  and  develop  a  method  to  pieilul  atleimaliou  a! 
any  fret|uenc\ .  given  common  sediment  properties. 


RI-.SIJ  L  IS 

In  si/ii  measurements  of  vhiiuI  velocity  and  attenuation  vvete  made  in  the  sea 
floor  oil  San  Diego  from  the  tesearch  submersible  lift  /»%/.// .  and  by  scuba  diving  in 
shallow  water,  during  the  period  |‘><.<i-7(l  Me.isuiementswerem.ule.il  '.and 
14  k  1 1  /  \dditional  measurements  are  reported  front  the  l‘b*'  hn\it  program  I  It  esc 

measiiiemenls.  and  others  trout  the  literature,  allowed  analysis  of  the  relationship  between 
attenuation  anti  Iretpiency.  and  study  ot  appropriate  visioelastu  models  vvliuh  van  be 
applied  to  saluraletl  seilnnents.  \n  understanding  ol  the  telalioiiships  between  attenua¬ 
tion  and  olliet  physiv.il  pioperlies.  such  as  yrtaiit  si/e  and  porosity,  allows  study  ot  the 
c.tuses  of  atleiiuaiion.  ii|x »n  which  predictions  van  be  based 

I  he  follow  my!  are  important  conclusions  cuiivcinmg  sound  attenuation  and  t elated 
properties  of  sea-floor  sediments. 

•  A  method  is  dented  which  allows  predivtion  ol  attenuation,  given  sediment 
mean  jtrain  si/e  or  |>orosity 

•  Atleimaliou  is  related,  approximately,  to  the  lint  |»*wei  ot  fictpicncy.  that  is. 
m  the  valuation,  a  A/",  n  is  about  I.  where  a  is  attenuation  of  coinpressioii.il 
t  sound  I  waves  in  decibels  per  metei.  A  is  a  voustaiit  which  applies  to  a  men 
sediment  type.  I.  the  Ireipiency  in  kilohertz,  and  u  is  ihe  exponent  of  fie<|uemy. 

•  I  here  is  no  significant  dependence  of  sound  velocity  on  licipicncy  t“  velocity 
dispersion”)  from  a  lew  her!/  to  the  megahertz  i.ingc. 


•  A  viscoelastic  model  is  lecommeiHleil  in  which  Ihe  I  ante  clast u  moduli,  p  and 
X.  aie  replaced  by  complex  moduli  l ft  ♦  /) /»  and  iX  ♦  /X'  I.  m  which  g.  X.  and 
density  govern  wave  velocity  and  Hi’  and  i\'  govern  energy  ila  in  pi  nr*  In  this 


in 


model,  alter  eliminating  negligible  factors.  energy  damping  (as  expressed  by  Ifie 
<| ik! lit >  factor.  Q.  the  specific  allemialion  factor.  I  '<J-  or  the  logarithmic  dccrc 
mentt  is  independent  of  frequency  (in  the  range  of  most  interest  in  underwater 
acoustics  t.  and  linear  attenuation  I  expressed,  lo*  example,  in  decibels  per  meter) 
is  proportional  to  the  first  power  of  frequency,  and  velocity  dispersion  is  negli¬ 
gible  or  absent. 

•  I  nergy  losses  due  to  intergrain  friction  appear  to  be.  by  far.  the  dominant 
cause  ol  attenuation  in  water-s.it ur.iled  sediments.  Viscous  losses  due  to 
relative  movement  of  |x»re  water  and  mineral  frame  are  piobably  negligible 


K  I- COM  M  r.NI>  AT  ION  S 

Additional  in  \iln  measurements  should  be  made  of  .inclination  m  sca-llooi  sedi¬ 
ments.  in  various  sedimentary  environments,  and  at  low  frequencies  t.i  lew  licit/  to 
5  kll/l.  Such  measurements  will  allow  tat  frequencies  of  most  interest t  1 1 1  better 
definition  of  the  relations  between  allemialion  and  frequency;  l2>  impiovcmcnt  of  the 
metliod.  promised  herein,  of  prediction  of  attenuation:  and  i.*i  determination  of  envi¬ 
ronmental  differences,  if  any.  of  attenuation  m  various  common  sediment  types. 

Ihe  gradients  of  attenuation  with  depth  m  thick  sediment  sections  and  the  filler- 
effects  of  various  sediment  layers  should  be  measured  and  studied. 
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Rl  I  I  RI  NCI  S 


INTRODUCTION 


In  the  sea  floor  off  San  Diego,  riming  the  period  l%(>  to  |l)7().  in  situ  measure¬ 
ments  were  marie  of  the  attenuation  of  compressional  (sound)  waves  in  marine  sediments. 
Measurements  in  deep  water  were  made  from  aboard  a  deep-diving  submersible.  In  shallow 
water,  scuba  divers  were  used.  In  this  program,  probes  were  inserted  into  the  sea  floor 
and  vompressional-wave  velocity  and  attenuation  (hereafter  called  “velocity”  anti  "attenua¬ 
tion")  were  measured  at  three  frequencies  (3.5.  7.  and  14  kl  1/ )  without  removing  the 
probes  from  the  sediment.  This  report  also  includes  a  few  measurements  at  25  kll/.  from 
the  bathyscaph  Trieste  program  in  l‘)62.  and  re-evaluated  measurements  of  attenuation  at 
100  kll/  from  an  earlier  scuba  diving  program  (Ref.  I  ). 

In  situ  measurements  from  the  above  programs  and  other  laboratory  and  in  situ 
measurements  of  attenuation  found  in  the  literature  allow  study  of  the  dependence  of 
attenuation  on  frequency  and  dispersion  (if  any)  of  compressional-wavc  velocity. 

The  relationships  between  frequency  and  compressional-wavc  velocity  and  attenua¬ 
tion  have  important  implications  in  forming  parameters  for  permissible  theoretical  elastic 
or  viscoelastic  models  for  water-saturated,  porous  sediments. 

The  objectives  of  this  study  are  to:  ( I )  report  in  situ  measured  values  of  velocity, 
attenuation,  and  associated  physical  properties  (/.<’..  density,  porosity,  grain  si/el  of 
marine  sediments;  (2)  assemble  and  analyze  pertinent  literature  data  on  attenuation; 

(3)  discuss  the  relationships  between  frequency,  velocity,  attenuation,  and  other  physical 
properties;  (4)  discuss  the  causes  of  attenuation  in  saturated  sediments;  (5)  discuss  elastic 
and  viscoelastic  models  which  can  be  applied  to  marine  sediments;  and  (6)  suggest  a  method 
to  predict  attenuation  in  marine  sediments,  given  frequency  and  common  physical  proper¬ 
ties.  such  as  porosity  and  grain  size. 


I  XPI  RIMI  NTAL  FINDINGS 


INTRODUCTION 

Three  sets  o  1/7/  situ  measurements  (made  by  the  writer)  of  the  attenuation  of 
compressional  waves  are  listed  in  fables  I  and  2  and  plotted  versus  frequency  in  f  igs.  I 
and  2.  In  addition,  some  in  situ  and  laboratory  measurements  from  the  literature  are 
listed  in  fables  3  and  4.  and  plotted  in  Fig.  2.  .1  or  comparison  and  later  discussions,  all 
values  are  listeil  (many  re-computed)  in  decibels  per  meter.  The  equations  showing  de¬ 
pendence  of  attenuation  on  frequency  are  listed  in  the  form 
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I  able  2 .  In  Sim  Attenuation.  Velmity. ami  Ollier  Physical  Properties of Sediments OH  Sail  Diego,  1056  I’r  igiani  * 


Sediment  t  ype 

No.  ol 
Stations 

1  Water 
Depth 
tint 

Moan  (train 

1  ti.iiin  tor 

i 

" 

Density 

(gin/em'l 

I'nrnsity 
'  1 1 

Velocity 

(in/soo) 

Altenuatiiin  idll/ml 
at  100  kll/ 

(,*> 

mill 

Sand 

t'oaiso 

i 

JO 

0.70 

0.578.1 

2.080 

.18..! 

1752 

5.1  1 

Medium 

i 

10 

1  00 

0.25 1 7 

2.000 

40.0 

16.10 

4  7  .1 

1  ino 

17 

III 

2.46 

0.1817 

1.026 

46.7 

1684 

52.1 

Very  lino 

4 

ll< 

•  lb 

0  1110 

1.0.18 

47.4 

1667 

55.0 

Sandy  nil 

2 

1.1 

4  2.1 

0.05.1.1 

1  860 

51.2 

1610 

74. .1 

Silly  sand 

17 

5.07 

0.0208 

1  680 

61  .1 

15.17 

60.0 

Sill 

1 

Is 

5.5.1 

0.0216 

1 .600 

60  0 

1465 

15  0 

Saml-vill'i  1  ii 

1 

17 

5.05 

0.0.102 

1.721 

58.0 

1400 

50.7 

Clayey  sill 

1 

*  22 

<>.10 

00146 

1  600 

65.6 

1464 

18.0 

'Averaged  values  lor  number  ol  nations  indicated  All  stations  occupied  l<y  setil'a  ilivers. 


where 


a  is  attenuation  of  compression;!  1  waves,  db/m 
A  is  a  constant 
./  is  fret)iiency,  kll/ 
n  is  the  exponent  of  frequency 


/.V  SITU  M  li ASUR  I:M  ENTS 
1966-1970  Program. 

The  equipment  am!  methods  used  to  make  in  situ  measurements  of  velocity  and 
attenuation  from  a  submersible  have  been  described  and  diagrammed  in  a  previous  report 
(Ref.  2).  In  general,  the  equipment  consisted  of  three  stainless-steel  probes.  7  cm  in 
diameter,  fastened  to  a  2-m-long,  rigid  beam  in  such  a  manner  that  when  the  beam  is  on 
the  sea  floor,  the  probes  are  inserted  a  variable,  pre-set  depth  into  the  sediment.  This 
depth  varied  from  30  to  60  cm  during  the  several  submersible  dives.  Three  barium- 
titanate  transducers  were  used,  one  as  a  sound  source  and  two  as  receivers.  Velocity  was 
determined  by  measuring  travel  time  over  a  I -in  path  between  the  receivers.  Attenuation 
was  measured  in  decibels  relative  to  that  in  the  bottom  water  (assumed  to  be  zero  for 
I  m).  Both  velocity  and  attenuation  were  measured  at  3.5.  7.  and  14  kll/  without  dis¬ 
turbing  the  probes.*  Coring  tubes  attached  to  each  end  of  the  rigid  beam  obtained  sedi¬ 
ment  sii tuples  for  laboratory  analyses.  Sediment  properties  were  obtained  by  standard 
laboratory  procedures.  Variations  and  inter-relationships  of  these  properties  as  well  as 
sediment  nomenclature  were  discussed  in  Kef.  3. 

*  I  ho  submersible  Ihi'i'Mur  4000,  from  which  mini'  ol  l  ho  measurements  wore  iii.iilo.  could  not  always  insert 
Iho  probes  to  ntasiniurn  poll ot ration  roipnroil  lo  obtain  measurements  at  .1.5  kll/. 


3 


Frequency,  kHz 


figure  1.  In  sill i  measurements  of  coiupressional-wa'c  ittcnuation  oft  San  Diego 
Symbols-  diamonds— 1966-1970.  Table  I:  inserted  trungk's-1962.  I.iblc  I: 
triangles— 1956,  Iablo  2.  loiters  indicate  sediment  type:  ( '-coarse  sand: 
M-modium  sand:  I  -fine  sand:  V-very  fine  sand:  SSI-sandy  silt:  SSA-siltv  sand: 
SSC-sand-silt-clay;  S-silt;  <  S-cla yes  silt,  line  labeled  "r"  indicates  slope  of 
any  line  having  a  dependence  of  attenuation  on  the  first  power  of  frequencs . 
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l'abli.'  3.  Laboratory  Attenuation  amt  Other  Physical  Properties 
of  Sediments  Off  San  Diego  I  rom  Kef.  IK.* 


ScJimcnt  type 

No.  of 

Mean  (train 

Diameter 

Density 

tgin/cni') 

Porosity 

Velocity 

23“C 

<ni/scc) 

Attenuation 

Stations 

0 

nun 

(  ) 

dB/m 

kit/ 

Sunil 

Medium 

13 

1.44 

0.36H6 

2.004 

39.1 

1737 

9.4 

26.6 

line 

13 

2.49 

0.1 7K0 

1 .9611 

43  6 

1693 

9.6 

26.(1 

Silly  valid 

1 

23/9 

0.1259 

1 .796 

53.2 

1551 

20.6 

23.9 

Sandy  sill 

4 

4.57 

0.0421 

1.737 

57.4 

1541 

11.8 

23.6 

Sill 

K 

5.27 

0.0259 

1.653 

615 

I54K 

7.8 

24.2 

Clayey  sill 

4 

6.31 

0.0126 

1  605 

65.5 

1508 

2.8 

23.2 

Clayey  sill 

1 1 

7.51 

0.0055 

1  336 

KO.O 

i 

1495 

16 

22.8 

*  A veraped  values  for  number  of  stations  indicated .  velocity  and  attenuation  lor  the  second  reso¬ 
nant  mode  (201 )  only. 


The  some  probes  were  used  in  shallow  water.  Operating  from  a  small  boat,  divers 
observed  the  sea  floor,  took  samples,  and  helped  insert  the  probes  into  the  sediment. 

The  results  of  the  above  measurements  and  subsequent  sediment  laboratory  analy¬ 
ses  are  listed  in  fable  I.  and  the  attenuation  versus  frequency  data  are  plotted  in  I'igs.  I 
and  2. 

1962  Program 

In  1 96  2.  velocity  and  attenuation  measurements  were  made  at  25  kHz  from  the 
bathyscaph  Trieste  and  from  a  small  boat  in  shallow  water.  The  equipment,  methods, 
and  results  of  the  velocity  measurements  were  reported  in  Kef.  4.  The  probes  used  for 
these  measurements  were  an  earlier  model  of  those  used  in  the  1966-70  program;  field 
methods  were  the  same. 

The  few  attenuation  measurements  made  during  this  program  (not  previously  re¬ 
ported  tare  listed  in  Table  I  anil  plotted  versus  frequency  in  figs.  I  and  2. 

1956  Program 

In  1956.  sound  velocity  and  attenuation  measurements  were  reported  for  a  num¬ 
ber  of///  situ  and  laboratory  measurements  in  shallow  water  off  San  Diego  (Kels.  I  and  5). 

In  situ  velocity  measurements  were  made  at  100  kllz.  The  equipment  and  methods 
have  been  reported  (Kef.  I  ).  In  general,  the  equipment  consisted  of  two  probes  secured 
through  a  rigid  base  so  that,  when  inserted  6  in.  (15  cm)  into  the  sea-floor  sediment  (by 
divers),  velocity  and  wave-amplitude  measurements  could  be  made  over  a  1  -ft  (50  cm) 
path  length.  At  the  same  stations,  small  containers  (10  cm  long  and  5  cm  in  diameter) 
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lahle  4.i  Attenuation  and  ( Mirer  Properties  of  Naluial.  Saturated  Sediments  tl  nun  the  I  iteuture)  * 


Sc  l  imcn  t  I  vpc 

Mea 

Dia 

i,v 

l  drain 

meter 

iiiiii 

Density 

(gm/citi'l 

Porosity 

1) 

\  ehvily 

25  C 
(m/sec) 

(t* 

»r‘ 

1  me  sand 

2  5 

(1  177 

<1  99) 

411) 

(1742) 

0  15 

1  26  •  0  14 

Veiy  fine  sand 

5  <15 

(1121 

(1  95) 

44  4 

(1716) 

(12  7 

1.17  •  0  14 

Very  line  \.md 

5.N 

(1072 

(1  97) 

42  5 

(1614) 

0  56 

1  (HI  •  0  01 

Silt 

5  42 

11.025 

1 1  901 

46  ti 

(1622) 

0  50 

1  05  •  o  15 

Medium  sand 

1  52 

(1.4(11 

1  99 

|40  H) 

11X15) 

0.164 

1 

Silty  wind 

1  <2 

1X2  2) 

11540) 

0  151 

1 

(  lay  ey  Mil 

1  20 

(X9.2) 

(155(1) 

0  1 IX 

1 

(  lay  -silt 

115 

192  1) 

(1525) 

0075 

1 

Silty  cl.iy 

114 

<95  5) 

11525) 

0072 

1 

•Mud" 

(7  5) 

<0.006) 

1  52 

76.0 

(14X0) 

0  1)66 

1 

(Clayey  silt) 

(7  <>l 

(I)  1)115) 

tKO.O) 

0  049 

1 

(Clayey  silt) 

0.066 

Sand -sill-clay 

(»  45 

(1.012 

1.75 

59.5 

1527 

(0  595) 

(0  99)  1 

Clay 

9  95 

0.001 

(1  451 

74.X 

( 1  (i  1 5 1 

51  9  dlt/m  at  46 X  kll/ 

54  dB/m  at  250  kll/ 

(Clayey  silt) 

(7.5  > 

(75.0) 

( 1405) 

26  2  dB/m  at  250  kll/ 

... 

(1621) 

0.2.  d B/iii  at  1  kll/ 

Null's 


I 

I 

I 

I 


5 

4a 

4t< 

5 

<> 

7  a 
7I> 
X 


*V  allies  in  parentheses  eomputeil  or  estinuited  In  the  writer  of  this  report 

1  In.  a  kj  ",  where  attenuation,  a.  is  in  illt/in:  frequency,/.  in  kll/;  k  is  a  eonstanl.  anil  n  is  t lie  expo¬ 
nent  of  Ireipieney. 

Notes: 


I .  Met  aim  ami  McCann  (Kef.  X);  in  situ,  land  and  heaeh;  velocity  given  at  20  1';  5  to  5(1  kll/. 

2  Mel  eroy  and  Del  oaclt  (Kef.  20); laboratory;  beach  sand  and  St.  Andrew  Hay.  I  lorida;  15  to 
1500  kll/ .  velocity  from  measured  ratio. 

.1.  Wood  and  Weston  (Kef.  9); in  situ,  tidal  mud  Hal.  I  nisworlh.  I  upland.  4  to  50  kll/;  velocity  from 
measured  ratio. 


4.  Hen  lie  It  (Kef.  I  2 1;  in  situ  i  reflection  teehimpie).  North  Atlantic  and  Mediterranean;  4a  12  kll/; 

4h  40  to  900  II/. 

5.  lewis  (Kef.  \A);insiln  (best  measurements),  shallow  water  oil  Puerto  Rico;  5  to  50  kll/. 

<>.  McCann  and  McCann  (Ref  X);  laboratory .  eight  core  samples.  North  Atlantic,  carbonate  >  5" ; 

46  X  kll/. 

7a.  Dlmiska  ( Kef.  10);  in  situ,  shallow  Baltic  Sea;  velocity  from  measured  ratio  (Station  52). 

7b.  I  lonska  (Kef.  10).  in  situ,  shallow  North  and  Baltic  Seas;  average  for  all  stations  with  velocity  ratios 
less  than  1 .1)1);  250  kll/.  velocity  from  measured  ratio. 

X.  Schirmcr  (Kef.  I  I ); in  situ;  shallow  Baltic  Sea  ( I  station),  velocity  from  measured  ratio. 


wore  inserted  by  h;tnd  into  the  sediment  lor  subsequent  laboratory  measurements  of 
velocity  and  attenuation  by  means  of  a  resonant-chamber  method  (Refs.  (i  and  7). 

In  Ref.  I.  velocities  were  reported  for  both  the  in  situ  probe  and  laboratory 
resonant-chamber  methods.  Values  of  attenuation  (at  23  to  40  kll/)  were  reported  for 
the  resonant-chamber  method. 

During  the  in  \ittt  probe  measurements  at  100  kll/..  the  amplitude  of  the  received 
wave  was  measured  on  the  oscilloscope  face  in  both  bottom  water  and  sediment;  three 
to  five  separate  measurements  were  made  at  each  station.  These  amplitudes  were  reported 
(Ref.  I.  Table  II)  as  “20  lop  (amplitude  water/amplitude  sediment)”,  and  not  as  linear 
attenuation,  because  of  uncertainties  concerning  impedance  loading  on  the  transducers  in 
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|jhlc4h  AttcmulKHi  in  Sediments  jnd  SediiticiilJiv  StulJ 
<1  i»ni  flic  I  ilctjluic)  * 


Sediment  1  v  pc 

A* 

«' 

1  ICt|UCIk  V 

III/) 

Nok*\ 

CUt  sand 

^  (MW  3 

1 

50  to  400 

1 

Malci-saluiatcd  clay 

(1  336 

1 

ISCIMDK  1 

2 

Scdinicntaty  strata 

II.J43 

1 

5  to  SO 

} 

Sedimentary  xiuii 

O.l  t** 

1 

IscIMIIk  ) 

4 

ScdimcnUiy  xtuu 

(MI35 

' 

(xcixinii ) 

5 

*<  Seismic  I:  assumed 

Mn.  <i  kfn.  win  tc  attenuation.  u.  is  in  UH/m:  lie- 


ipicnc)  .  hi  kll/;  k  is  a  constant.  ami  n  is  the  exponent  ol  Ire 
l|IICIIC) 

Notes: 

1 .  I  ullus  and  Keid  (Kid  I  3).(iiilf  td  Mexico  coaxial 
sediments; depth  internal  3  lo  34  in. 

2.  /Iiadin  (in  Kcl  1 5);  depth  interval  0  to  JO  in. 

3.  /enixtov  (Kef.  16). 

4. 1  pinatyeva  i  t  til  (in  Kef.  16) 

5.  Her/on  (in  Kef.  16). 


the  two  media.  Re-evaluation  of  the  method  and  results  indicates  these  data  van  he 
identified  as  attenuation  in  the  sediments.  The  basic  data  have  been  re-analy/ed  and 
extrapolated  to  I  in  and  are  listed  in  the  form  of  decibels  per  meter  in  Table  2:  attenua¬ 
tion  versus  frequency  is  plotted  in  Figs.  1  and  2. 


MEASUREMENTS  FROM  THE  LITERATURE  UN  SITU) 

In  situ  measurements  of  sound  attenuation  in  marine  sediments  are  relatively  rare 
compared  with  measurements  of  other  properties  of  sediments.  Those  of  the  writer  are 
noted  above.  In  this  section,  brief  resumes  of  other  in  situ  measurements  in  the  literature 
will  be  presented.  Included  are  in  situ  measurements  on  land  in  some  water-saturated 
sediments  and  sedimentary  rocks,  which  may  be  pertinent  to  lower  layers  of  the  sea  floor 
or  to  later  discussions.  These  attenuation  measurements  are  plotted  versus  frequency  in 
Fig.  2,  and  listed  in  Tables  3  and  4a,  b.  For  those  studies  in  which  equations  are  given 
Refs.  X  and  ()).  points  are  plotted  in  Fig.  2  for  the  lowest  and  highest  frequencies 
and  an  intermediate  frequency. 

Wood  and  Weston  (Ref.  ‘))  made  in  situ  measurements  of  attenuation  (4  to 
50  kll/)  in  tidal  mud  flats  in  Hmsworth  Harbor.  England.  Maximum  distance  between 
probes  was  1X0  ft  (55  m):  probe  penetration  was  about  0.5  m. 

(.Honshu  (Ref.  10)  and  Schirmer  (Ref.  1  1  )  measured,  in  situ,  compressional  velocity 
and  attenuation  in  shallow-water  sediments  in  the  Baltic  Sea  at  sediment  depths  of  about 
2  m.  Ulonska’s  measurements  were  over  a  0.5-m  path  at  250  kll/:  Schirmcr's  over  a  I  2-m 
path  at  about  I  kll/. 


X 


McCann  and  McCann  (Ref.  8)  reported  in  situ  measurements  in  water-saturated 
land  and  beach  sediments  at  a  mean  sediment  depth  of  2  m.  Frequencies  were  between 
5  and  50  kHz. 

Bennett  (Ref.  12).  using  reflection  techniques  from  a  surface  vessel,  computed 
values  of  attenuation  at  12  kHz,  and  between  40  and  ()()0  Hz,  for  sediments  in  the  North 
Atlantic  and  Mediterranean. 

Tullos  and  Reid  (Ref.  13)  measured  attenuation  in  Ciuifof  Mexico  coastal  sedi¬ 
ments  (clay-sand)  of  Pleistocene  age  in  three  depth  intervals  to  1000  ft  (305  m);  fre¬ 
quencies  were  between  50  and  400  Hz.  Measurements  in  the  depth  interval  of  7.5  to 
1 10  ft  (2  to  34  m)  are  listed  in  Table  4b. 

Lewis  (Ref.  14)  measured  attenuation,  in  situ,  at  one  station  in  shallow  water  off 
Puerto  Rico  by  inserting  probes  into  the  sediment:  frequencies  were  5  to  50  kHz.  The 
data  were  scattered,  but  Lewis  believed  they  represented  a  dependency  of  attenuation  on 
frequency  to  the  first  power. 

There  are  several  studies  from  Russian  literature  listed  in  Table  4b.  Zliadin  (in 
Ref.  15)  reported  attenuation  in  “water-saturated  clay”  in  the  depth  range  of  0-20  m. 

No  frequency  is  given,  but  is  presumed  to  be  in  the  seismic  range.  Zemtsov  (Ref.  16) 
reports  his  own  studies  of  attenuation  in  the  frequency  range  of  5  to  50  Hz  in  “sedi¬ 
mentary  strata”  and  notes  attenuation  versus  frequency  equations  from  Fpinatyeva  and 
Berzon  in  similar  materials;  values  are  shown  at  30  Hz.  in  big.  2. 


MEASUREMENTS  FROM  THE  LITERATURE  (LABORATORY) 

It  is  possible,  but  difficult,  to  make  valid  laboratory  measurements  of  attenuation 
in  natural  or  artificial  water-saturated  sediments.  Even  when  valid  measurements  are 
made  in  some  experimental  materials,  the  results  cannot  always  be  extrapolated  to  natural 
marine  sediments.  Some  of  these  difficulties  are  noted  below. 

One  of  the  chief  difficulties  in  measuring  attenuation  in  “water-saturated”  labo¬ 
ratory  samples  is  air  entrapped  in  the  pore  spaces  of  artificial  sediments,  or  air  or  gas 
from  organic  materials  in  natural  sediment  samples  transported  to  the  laboratory.  When 
air  or  gas  is  present  in  “water-saturated” sediments,  the  measured  attenuation  is  apt  to  be 
too  high  and  velocity  too  low  (depending  on  frequency  and  bubble  size).  The  experi¬ 
ences  of  Wood  and  Weston  (Ref.  l))  with  sands  and  muds  are  instructive  in  this  matter. 

Some  artificial  sediments  are  not  truly  analogous  to  natural  sediments.  Lor  ex¬ 
ample,  kaolinite  in  distilled  water  containing  a  deflocculent  acts  as  a  Newtonian  fluid 
without  rigidity,  whereas  in  a  flocculated  structure,  rigidity  is  present  (Ref.  17).  As  dis¬ 
cussed  in  a  later  section,  the  relation  of  attenuation  with  frequency  differs  in  the  two 
cases.  A  favored  laboratory  material  is  round -grained  quartz  sand  (e.g..  St.  Peter’s  sand). 


This  material  lias  less  rigidity  and  attenuation  than  the  more  common.  angular-grained 
natural  sands  (discussed  below I. 

Some  laboratory  measurements  in  natural  marine  sediments  of  parlieular  sigmli- 
eanee  to  the  present  study  are  noted  and  briefly  discussed  in  this  section. 

Shumwuy  (Kef.  IS)  used  a  resonant -chamber  technique  to  measure  attenuation 
and  velocity  at  various  frequencies  between  20  and  40  kll/.  The  basic  theory  and  equip¬ 
ment  were  reported  in  Refs.  (>  and  7.  These  measurements  have  been  widely  referenced, 
and  values  extrapolated  to  other  frequencies. 

Unfortunately.  Shumway's  measurements,  alone,  cannot  be  used  to  determine 
frequency  and  attenuation  relationships,  although  in  a  number  of  samples,  measurements 
were  made  at  two  or  three  frequencies.  A  detailed  examination  of  the  published  data 
and  the  unpublished  records  of  measurements  (in  the  writer's  files)  reveals  numerous 
experimental  errors.  These  errors  are  evident  for  the  following  reasons:  ( I )  in  individual 
samples  where  attenuation  was  measured  at  two  or  three  frequencies,  there  are  In  eases 
(in  Kef.  IS)  in  which  attenuation  dvitcascil  with  increasing  frequencies,  a  result  incon¬ 
sonant  with  all  theory  and  other  experiments  in  similar  materials,  and  (2)  even  with  these 
In  cases  eliminated,  computations  of  the  frequency  dependence  of  attenuation  (in  the 
form:  a  =  kf")  revealed  that  the  exponent  ii  varied  in  similar  sediment  types  (having 
similar  physical  properties)  as  follows:  sand  <0.0  to  3.1 ).  sandy  silt  and  silly  sand  (O.S 
to  2.0).  and  clayey  silt  (0.0  to  3.4).  Neither  theory  nor  other  experiments  in  similar 
materials  support  these  variations. 

The  causes  of  the  above-noted  experimental  errors  are  unknown.  Shumway's 
samples  were  excellent.  Most  were  taken  in  sirn  by  divers  and  never  exposed  to  air  prior 
to  measurement.  These  errors  may  be.  in  part,  related  to  the  method.  Shumway  noted 
(Kef.  0.  p.  318)  that  the  resonant-chamber  method  assumes  samples  with  no.  or  negligible, 
rigidity.  This  may  be  true  in  some  shallow-water  silt-clays,  but  not  in  sands  and  silty 
sands.  (See.  for  example,  data  in  Ref.  10.)  In  addition,  not  all  resonant  modes  were  of 
equal  reliability.  Toulis  (Kef.  7)  suggested  that  the  second  resonant  mode  (201  )  was  most 
reliable. 


However.  Shumway's  attenuation  data  (Ref.  18)  can  probably  he  used,  selectively, 
by  averaging  values  for  general  sediment  types;  hopefully,  experimental  errors  will  cancel 
out.  Tor  comparison  with  the  in  siln  measurements  off  San  Diego  (presented  herein)  and 
with  other  values  from  the  literature,  the  writer  averaged  attenuation  values  (and  other 
properties)  at  the  second  resonant  mode  for  various  sediment  types  off  San  Diego;  these 
are  listed  in  Table  3  and  plotted  in  Tig.  2. 

McLeroy  and  DeLoach  (Ref.  20)  measured  (in  the  laboratory)  sound  speed  and 
attenuation,  from  15  to  1500  kll/.,  in  natural  sediments  from  sites  in  St.  Andrew  Day. 
Panama  City.  Florida,  and  in  sand  from  a  (iulf  of  Mexico  beach.  They  noted  no  indica¬ 
tion  of  entrapped  gases.  These  measurements  are  listed  in  Table  4a  and  plotted  in  l  ig.  2. 
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I  lie  unusually  low  densities  :itul  computed  (by  the  writer)  high  porosities  of  the  Hoy  .edi- 
inents  lire  not  typical  of  most  open-ocean  sediments, 

I  he  excellent  study  of  attenuation  in  marine  sediments  hy  McCann  and  McCann 
(Kef.  K)  lias  been  noted  in  the  section  on  in  \iln  measurements.  These  writers  also  mea¬ 
sured  attenuation  at  4oX  k 1 1 /  in  sediment  cores  from  the  North  Atlantic.  Most  ol  these 
sediments  contained  appreciable  amounts  of  calcium  carbonate.  There  is  some  indication 
that  attenuation  in  this  type  of  sediment  may  be  higher  than  in  non-calcarcous  day-silts 
with  similar  grain  si/cs  and  porosities,  because  calcareous  sediments  are  not  included  in 
this  report,  an  average  value  of  attenuation  for  eight  clay  samples  with  less  than  5  percent 
carbonate  is  listed  in  Table  4a  and  plotted  in  Tig.  2. 


DISCUSSION  AND  CONCLUSIONS 


INTRODUCTION 

The  study  of  elastic  and  viscoelastic  models  which  can  be  applied  to  dry  and 
saturated  porous  rock  and  sediments  has  lonccrncd  a  large  number  of  scientists  and 
engineers,  Many  of  these  studies  have  considered  the  extent  to  which  these  media  can 
be  described  by  the  equations  of  llookean  elasticity  or  by  those  of  the  Kelvin-Voigt. 
Maxwell,  or  other  viscoelastic  models.  It  is  somewhat  surprising,  therefore,  to  discover 
such  a  wide  diversity  of  models,  equations,  and  opinions  on  such  an  important  subject. 

In  the  case  of  water-saturated  natural  sediments,  the  reason  for  this  diversity  appears  to 
be  that,  until  leeently.  experimental  evidence  capable  of  restricting  model  parameters  has 
been  scarce.  In  the  following  discussion  (unless  otherwise  noted)  it  is  assumed  that  the 
medium  is  a  porous,  uncemented  mineral  structure,  fully  saturated  with  water.  The  stress 
is  that  of  a  eompression.il  or  shear  wave  of  low  amplitude,  and  strains  are  of  the  order  of 
I  O'6,  or  less.  Wavelengths  are  much  greater  than  grain  si/e.  otherwise  Rayleigh  scattering 
can  occur:  and  attenuation  is  related  to  the  fourth  power  of  frequency  (see.  for  example. 
Ref.  :i>. 

The  equations  of  llookean  elasticity  do  not  account  for  energy  damping.  Conse¬ 
quently.  an  adequate  model  must  be  anelastic  if  energy  damping  is  considered.  In  the 
selection  of  an  appropriate  anelastic  model,  a  critical  factor  is  the  extent  of  relative  move¬ 
ment  of  pore  water  and  mineral  particles.  If  the  pore  water  moves  significantly  relative 
to  mineral  structure,  then  viscous  damping  and  velocity  dispersion  must  be  considered. 

If  the  pore  water  does  not  move  significantly  with  respect  to  the  solids,  there  is  negligible, 
or  no.  velocity  dispersion,  and  energy  damping  is  not  dependent  on  viscosity  of  pore 
water  anil  permeability  of  the  mineral  structure.  The  dependence  of  energy  damping  on 
frequency  is  different  in  the  two  cases.  Consequently,  two  critical  parameters  for  anelas- 
lic  models  are  velocity  dispersion,  if  any.  and  Ihe  dependence  of  energy  damping  on 
frequency.  In  the  next  section,  the  experimental  evidence  concerning  these  two  factors 
will  be  reviewed.  In  following  sections,  a  particular  anelastic  model  is  recommended,  and 
the  relationships  between  attenuation  and  other  physical  properties  and  Ihe  causes  of 
attenuation  will  be  discussed.  In  the  last  section,  a  method  for  predicting  attenuation 
(given  grain  si/e  and  sediment  porosity)  will  be  discussed. 


SOME  PARAMETERS  Of  ELASTIC  AND  VISCOELASTIC 
MODELS  FOR  SATURATED  SEDIMENTS 


Introduction 

Recent  experimental  studies  in  wave  velocities  and  attenuation  and  in  dynamic 
rigidity  have  placed  important  restrictions  on  probable  (and  possible)  elastic  and  visco¬ 
elastic  models  for  water-saturated  sediments.  These  parameters  merit  a  more  extended 
review  because  of  their  importance  in  geophysical  studies.  Additionally,  the  concepts 
and  statements  of  this  paper  require  documentation.  These  restrictions  apply  to  porous 
sediments  saturated  with  water,  and  without  a  gas  phase,  when  wavelengths  are  much 
greater  than  grain  sizes.  The  frequency  range  is  from  a  few  hertz  to  several  hundred  kilo¬ 
hertz,  or  into  the  megahertz,  range.  These  restrictions  are  concerned  with  the  question  of 
the  dependence,  if  any,  of  velocity  on  frequency  (“velocity  dispersion"),  and  the  depen¬ 
dence  of  wave-energy  damping,  or  attenuation,  on  frequency.  In  the  following  sections, 
these  restrictions  will  be  documented  and  discussed. 


Velocity  Dispersion 

A  number  of  investigators  have  measured  compressional-  and  shear-wave  velocities 
in  rocks  (laboratory  and  in  situ).  Most  have  concluded  that  there  is  no  (or  negligible) 
measurable  velocity  dispersion  in  the  range  from  seismic  frequencies  into  the  megahertz 
range.  Examples  include  work  and  reviews  by  Wylie  ct  al.  (Ref.  22).  Birch  (Ref.  23). 
Peselnick  and  Outerbridge  (Ref.  24).  White  (Ref.  25).  and  Press  (Ref.  2b). 

Most  viscoelastic  models  requiring  a  dependence  of  velocity  on  frequencx  include 
movement  of  pore  water  relative  to  the  mineral  frame,  which,  in  turn,  causes  water  vis¬ 
cosity  and  sediment  permeability  to  be  important  factors.  Sands  have  relatively  high 
permeability,  large  grain  size,  and  interconnecting  pores:  consequently,  these  models  indi¬ 
cate  maximum  velocity  dispersion  in  sands.  It  is  therefore  instructive  to  examine  the 
experimental  evidence  regarding  velocity  dispersion  in  sands  across  wide  frequency  ranges. 

Because  permeability  is  a  factor  in  some  models,  it  is  pertinent  that  Wylie  ct  al. 
(Ref.  22 1  found  that  velocity  through  water-saturated  glass  beads  of  various  sizes  was  un¬ 
affected  as  permeability  varied  by  a  factor  of  4.6  X  I  O'4. 

No  velocity  dispersion  was  measured  in  clean,  round-grained  sands  in  the  laboratory 
by  Hunter  ct  al.  (Ref.  27)  over  the  range  7  to  73  kHz;  Hardin  and  Richart  (Ref.  28).  0.2 
to  2.5  Ml/;  Nolle  ct  al.  (Kef.  26).  200  to  1000  kll/;  and  Schbn  (Ref.  30).  20  to  64  kHz. 

In  round -grained,  pure  quartz  sands  in  distilled  water  at  porosities  of  36  percent, 
Shumway  (Ref.  18,  p.463)  and  Nolle  ct  al.  (Ref.  26)  measured  compressional  velocities 
of  1744  m/sec  (26  kll/)  and  1740  m/sec  (400  to  1000  kll/).  respectively. 


In  soil  mechanics  investigations  in  situ  in  sands.  !ow-l'rci|iioncy  vibrations  were 
used  in  studies  by  Barkan  (Ref.  31 )  and  Jones  (Ref.  32)  to  measure  shear-wave  velocities 
(Jones’  measurements  also  included  clay-silt).  No  velocity  dispersion  was  measured  in 
the  frequency  range  from  10  to  400  II/. 

Ideally,  to  test  for  velocity  dispersion,  measurements  should  be  made  on  the  same 
sample  in  the  laboratory,  or  in  situ  in  the  same  sediment,  by  merely  changing  the  fre¬ 
quency.  In  natural  sands  this  has  been  done  in  the  laboratory  and  in  situ.  In  the  labo¬ 
ratory.  McLeroy  and  DeLoach  (Ref.  20)  measured  a  ratio  (velocity  in  sediment/velocity 
in  water)  of  1.189  in  medium  sand  over  a  frequency  range  of  15  to  1500  kll/.  In  situ 
McCann  and  McCann  (Ref.  8)  reported  no  velocity  dispersion  in  fine  or  very  fine  sands 
between  5  and  50  kll/.  During  the  1970  measurements  off  San  Diego  reported  herein, 
special  tests  were  made  in  fine  sand  at  four  stations  to  determine  the  presence  or  absence 
of  velocity  dispersion  between  3.5  and  14  kll/.  The  results  indicated  no  measurable 
velocity  dispersion  between  these  frequencies. 

In  natural  saturated  sands  it  is  difficult  to  compare  velocity  measurements  in  dif¬ 
ferent  samples  or  from  different  stations  because  of  the  effects  on  velocity  of  grain  si/e 
and  shape,  porosity,  and  other  factors,  d  or  a  discussion  of  these,  see  Ref.  3.)  However, 
the  following  in  \itu  measurements  (corrected  to  23 V)  indicate  no  significant  velocity 
dispersion  in  natural  fine  sand  over  a  frequency  range  from  14  to  100  kll/: 
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In  very  fine  sand.  McCann  and  McCann  measured  a  velocity  of  I7I<>  (poiosity.  43.3  per¬ 
cent)  at  frequencies  between  5  and  50  kll/. 

In  high-porosity  silt-clays  in  the  laboratory  and  in  the  field,  eompressional  velocity 
is  usually  less  in  the  sediment  than  it  is  in  the  water,  and  no  velocity  dispersion  is  indi¬ 
cated  over  very  wide  frequency  ranges.  The  evidence  is  summarized  in  Table  5.  Velocity 
data  are  presented  as  the  ratio:  velocity  in  sediment/velocity  in  water.  As  discussed  in 
Ref.  3.  this  ratio  remains  the  same  in  both  the  laboratory  and  the  surficial  sediments  of 
the  sea  floor. 

Some  recent  experiments  in  artificial  clays  in  the  laboratory  have  implications  for 
studies  of  possible  shear-wave  velocity  dispersion.  Cohen  (Ref.  17)  measured  complex 
rigidity  (i u  +  in')  in  flocculated  kaolinite  in  distilled  water.  Both  g  and  in'  were  inde¬ 
pendent  of  frequency  in  the  range  8.6  to  43.2  kll/.  Hardin  and  Black  (Ref.  40).  using 
a  vibration  technique  to  measure  dynamic  rigidity  in  kaolinite  in  distilled  water,  demon¬ 
strated  that  dynamic  rigidity  was  independent  of  wave  amplitude,  in  their  samples,  at 
amplitudes  less  than  1 0'4,  and  independent  of  frequency  between  200  and  300  Hz.  No 
dispersion  in  dynamic  rigidity,  p,  indicates  no  dispersion  in  shear-wave  velocity. 
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I  attic  5.  Summary  of  Ratios  of  Compressional-Wavc  Velocities  ( I’j,  Sediment/ 1  ^  Water) 
in  lligh-l’orosity  Silt-Clays  at  Various  frequencies. 


Material 

Velocity  Ratio 

frequency 

Notes 

Kaolinite  in  distilled  water 

<1.97  to  0.99 

1  Mil/ 

1 

Kaolinite  in  distilled  water 

0.97  to  0.99 

9  to  42  kll/ 

s 

Deep-sea  elay  slurry 

0.97 

28  kll/ 

.2 

Sand-silt-clay 

0.98 

14  kll/ 

4 

Clayey  silt,  San  Diego  1  rough 

0.98 

22  to  40  kll/ 

.2 

Clayey  silt,  San  Diego  Trough 

0.98 

25  kll/ 

5 

Clayey  silt,  San  Diego  Trough 

0.98 

14  kll/ 

4 

Clay-silt.  St.  Andrew  Hay 

0.997 

15  to  1500  kll/ 

6 

Silt,  clayey  silt,  continental  shelf 

0992 

2  Mil/ 

7 

Silt-clays,  continental  shelf 

0  994 

2  Mil/ 

8 

Silty  elay.  continental  shell 

0.994 

200  kll/ 

9 

Deep-sea  silt-clay 

0.977 

2(1  to  200  11/ 

10 

Deep-sea  sill-clay 

0.980 

20  to  20011/ 

1  1 

Deep-sea  silt-clay 

0.985 

201)  kll/ 

9 

Deep-sea  elay 

0.975 

200  kll/ 

9 

Deep-sea  silt-clay  (1  micron) 

0.977 

400  kll/ 

12 

Deep-sea  silt-clay  (2  microns) 

0.987 

400  kll/ 

12 

Deep-sea  silt-clay 

0.986 

40(1  kll/ 

12 

"1  ow -velocity  layer.”  continental  shelf 

0.980 

250  kll/ 

14 

Notes: 

1 .  Crick  (Ref. 22);  ratio  dependent  on  concentration  of  solids. 

2.  Cohen  (Ref.  1 7);  ratio  dependent  on  concentration  of  solids. 

2.  Shu nns a)  (Ref.  18);  resonant  chamber;  laboratory 

4.  Hamilton,  this  report;  w  silu  measurements. 

5.  Hamilton  (Ref.  2);w  situ  measurements. 

6.  Mel  croy  and  Deloach  (Ref.  20); laboratory  samples  from  Si  Andrew  Hay.  I  la 

7.  Hieda  (Ref.  24).  lops  of  5  cores  off  Southern  California. 

H.  1  asswell  (Ref.  25);9  cores;  tidal  mud  Hat.  Southern  California. 

9.  Hamilton  (Ref.  2);  laboratory :  samples  from  North  Pacific. 

10.  I  ry  and  Raitl  (Ref.  26);  reflection  technique ;  deep  Pacific. 

1 1 .  limit/  and  l.wing  (Ref.  27);  relleetion  teehnnpie.  deep  Atlantic 

12.  Horn  cl  al  (Ref.  28);  laboratory; core  samples  from  North  Pacific 
1 2.  Schreiber  (Ref.  29);  tops  of  1(1  cores  off  Hawaii. 

14  Clonska  (Ref.  ID);  tu  silu  measurements;  shallots  North  and  Baltic  Seas. 


Hiimpton  (Ref.  41)  measured  compression;! I  velocity  and  attenuation  in  artificial 
sediments  and  reported  a  marked  velocity  dispersion  (4  to  (>  percent)  in  silt-clays  between 
frequencies  of  3  to  200  kHz.  In  his  low-frequency  measurements  (3  to  30  kll/)  in  floc¬ 
culated  kaolinite  in  distilled  water,  he  reported  (Ref.  41,  Tig.  I  I  )  anomalously  low 
velocity  ratios  (0.03  to  0.94)  and  unusually  high  attenuation  values  (Ref.  41.  p.  88(i).  The 
experimental  evidence  of  other  investigators  (e.g. ,  Refs.  17  and  33)  in  the  same,  or  similar, 
materials  and  at  the  same  frequencies,  in  both  the  laboratory  and  field  (see  Table  5).  does 
not  support  Hampton’s  reported  low  velocity  ratios  or  velocity  dispersion.  There  is  a 
possibility  that  his  day-water,  artificial  sediments  were  not  air  or  gas  free. 

In  summary,  the  following  can  be  concluded  in  regard  to  velocity  dispersion  in 
saturated  sediments,  in  sands,  a  number  of  studies  (cited  above)  have  reported  no  dis¬ 
persion  over  restricted  frequency  ranges.  These  studies  and  comparison  of  the  values  of 
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velocity  in  similar  sands  a  I  different  frequencies,  as  above,  indicate  that  velocity  disper¬ 
sion,  if  present,  is  negligibly  small  from  a  few  kilohertz  to  the  megahertz  range.  The 
evidence  (as  in  Table  5)  indicates  that  velocity  dispersion  ir  higher  porosity  silt-clays,  it 
present,  is  negligible  over  a  frequency  range  from  less  than  I  kHz  to  2  Mllz.  However, 
it  must  be  stated  (as  a  reviewer  pointed  out)  that  most  of  the  quoted  tests  of  velocity 
dispersion  were  made  over  only  an  order  of  magnitude  of  frequencies  or  less,  which  is 
not  necessarily  enough  to  show  dispersion.  In  other  words,  it  cannot  be  stated  on  the 
basis  of  present  experimental  evidence  that  velocity  dispersion  is  non-existent,  especially 
over  very  wide  frequency  ranges  (from  a  few  hertz  to  several  megahertz). 

Energy  Damping 

The  relationship  between  frequency  and  wave-energy  loss,  or  damping,  is  a  critical 
parinv’tor  in  nekvlHw  uf  m  appropriate  Hivelurtiv  uw4el  tw  any  moUimu  Rnwirt  review* 
have  summarized  a  large  number  of  laboratory  and  field  studies  of  wave-energy  losses  in 
rocks  in  which  the  specific  attenuation  factor.  I  !Q.  and  the  logarithmic  decrement  have 
been  shown  to  be  approximately  independent  of  frequency  over  a  range  of  at  least  10s  II/. 
(Refs.  25  and  42-45).  Attwell  and  Rainana  (Ref.  45)  included  some  sediment  data. 

Evidence  that  the  specific  attenuation  factor,  \/Q,  is  independent  of  frequency 
implies  that  attenuation,  a.  in  decihels/uuit  length,  increased  linearly  with  frequency. .) 
(t’.g..  Ref.  25.  p.  ()X).  Recent  summaries  of  work  in  this  field,  in  the  ease  of  compres- 
sional  waves  (Refs.  25.  43.  and  45).  indicate  that,  for  most  rocks,  there  is  a  small  varia¬ 
tion  around  {meant y  w  the  range  frequencies tt?  mtcrc-st  m  ttndvfwaVci  acoustic; 
and  marine  geophysics:  that  is.  in  the  relationship  a  =  A./'",  the  exponent  //  is  approxi¬ 
mately  I.  These  studies  included  dry  rocks  (usual  in  the  laboratory)  and  in  situ  measure¬ 
ments  in  rock  strata  which,  below  groundwater  levels,  are  saturated. 

Some  recent  studies  of  attenuation  in  saturated  sediments  which  are  pertinent  to 
the  present  report  have  been  noted  in  a  previous  section  and  their  results  presented  in 
tables  and  figures.  Of  special  interest  in  these  studies,  and  in  the  measurements  of  this 
report,  is  the  dependence  of  attenuation  on  frequency.  In  the  following  paragraphs  the 
current  experimental  evidence  on  this  subject  will  be  reviewed. 

In  examining  the  dependence  of  attenuation  on  frequency,  it  is  important  to 
recognize  that,  in  natural  marine  sediments,  each  set  of  measurements  and  resulting 
equation  (e.g. .  a  =  kf ")  is  apt  to  be  unique.  The  reasons  for  this  phenomenon  will  be 
discussed  below;  but  in  general,  the  important  variables  involve  the  varying  sediment 
structures,  porosity,  grain  size  and  shape,  inlerparticle  contacts  and  surface  areas  in  sands 
and  coarse  silts,  and  physicochemical  forces  (cohesion)  in  the  fine  silts  and  clays.  These 
variable  factors  result  in  much  scatter  of  measured  values  of  attenuation  in  similar  sedi¬ 
ments  (Tables  1-4),  and  consequent  scatter  in  computed  values  of  the  constant  A.  Also, 
because  of  the  difficulties  of  making  accurate  measurements  of  attenuation  in  natural 
sediments  (in  situ  or  laboratory),  it  is  to  be  expected  that  computed  values  of  the  expo¬ 
nent  of  frequency,  n,  will  vary.  As  a  result,  generalized  statements  about  the  exact  de¬ 
pendence  of  attenuation  on  frequency  must  be  qualified  and,  as  far  as  justifiable,  a 
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statistical  approach  followed  in  studying  data  (which  was  the  approach  of  Attwell  and 
Kamana,  Ref.  45). 

The  important  factor  in  predicting  attenuation  and  extrapolating  to  various  fre¬ 
quencies  is  the  exponent  of  frequency.  //  (in  the  equation  above).  Also,  this  exponent 
is  a  critical  parameter  for  selecting  appropriate  anelastic  models  which  can  he  applied  to 
saturated  sediments.  In  Tables  I  and  4a  are  values  of  //  for  measurements  made  at  two 
or  more  frequencies  in  a  single  sample  or  sediment  type  by  an  individual  experimenter. 

In  Table  4b,  values  are  given  for  some  low-frequency  measurements  in  sediments  or 
sedimentary  rock  sections  (some  of  the  reports  did  not  specify  the  degree  of  lithifica- 
tion).  Some  of  the  measurements  in  Table  4b  may  be  pertinent  to  lit hified  or  semi- 
lithified.  deeper  layers  in  the  sea  floor.  As  can  be  seen  in  Tables  I  and  4,  measured 
values  of  //  in  these  natural  sediments,  from  widely  scattered  geographic  areas,  vary 
closely  around  I  for  a  wide  variety  of  deepwater  and  shallow-water  sediments  and  over 
a  wide  frequency  range. 

Computations  were  made  (Table  6,  Part  A)  which  interrelated  /'//  situ  measure¬ 
ments  in  similar  sediments  from  various  experiments  conducted  by  the  writer  in  the  sea 
floor  off  San  Diego  (Tables  1.  2).  T his  allowed  extension  of  the  frequency  range:  fre¬ 
quencies  usually  included  3.5,  7,  14,  25,  and  100  kHz.  Off  Mission  Beach,  attenuation 
data  for  two  stations  came  from  Ref.  I.  Table  II.  The  in  situ  data  for  “Same  area  off 
Mission  Beach"  and  "San  Diego  Trough”  involved  only  two  or  three  stations  in  each 
area.  In  these  locations  the  values  of  n  (rounded  to  two  decimals  in  the  table*)  are  1 .007. 
1.024,  0.947,  and  0.9b9.  The  data  for  "General  sediment  type"  for  fine  sand  involved 
24  stations:  the  value  of  n  is  1.007. 

In  Table  6,  Part  B.  the  writer's  in  situ  measurements  (Tables  I.  2.  and  h.  Part  A) 
were  combined  with  averaged  laboratory  values  for  general  sediment  types  off  San  Diego 
(Table  3)  from  Ref.  IS.  The  value  of//  for  37  stations  in  fine  sand  is  0.992:  for  medium 
sand  at  lb  stations.  //  is  0.‘)S2. 

Because  of  the  larger  number  of  stations,  the  data  for  “General  sediment  type” 
in  Table  b  are  considered  more  reliable:  1.007.  0.992.  0.082.  As  in  previous  tables,  the 
values  of//  vary  closely  around  I,  with  the  more  reliable  data  at.  or  slightly  below.  1. 

In  reconciling  acoustic  theory  with  experimental  acoustic  energy  losses  of  sound 
incident  on  the  sea  floor,  it  is  usually  necessary  to  assume  values  of  attenuation  at  fre¬ 
quencies  of  interest  for  various  sediment  types  and  layers.  Such  reconciliation  has  been 
successful  in  several  studies  (frequencies  from  0.1  to  4  kHz)  in  which  a  first-power  de¬ 
pendency  of  attenuation  on  frequency  was  assumed  (c,e  .  Refs.  4b-49).  Cole  concluded 
that  the  first-power  dependency  can  be  extended  over  the  range  100  to  TOO  11/ 

Compressional  velocity  and  attenuation  data  from  Refs.  10  and  II  are  possibly 
significant  in  determining  the  relationship  between  attenuation  and  frequency.  Schirmer, 
at  a  water  depth  of  20  m  in  the  Baltic  Sea,  measured  attenuations  of  0.2  to  0.4  dB/m  at 
1.0  to  1.4  kHz;  an  average  velocity  ratio  (in  sediment/in  water)  of  1 ,08b  (1.05  to  1.14) 
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was  measured  in  a  eore  from  the  site.  Ulonska  (Kef.  10).  at  a  water  depth  of  22  m  in  the 
Baltic  Sea.  measured  (Station  32)  a  velocity  ratio  of  1.082  and  attenuations  of  50  to  58 
dB/m  at  250  kHz.  Both  sets  of  measurements  were  at  about  2  m  depth  in  the  sea  floor. 
The  velocity  ratios  indicate  the  sediment  types  were  about  the  same.  If  a  first-power 
dependency  of  attenuation  on  frequency  is  assumed,  extrapolation  of  Ulonska 's  average 
data  (54  dB/m  at  250  kHz)  to  1  kHz  yields  an  attenuation  of  0.22  dB/m  (Schirmcr 
measured  0.2  dB/m  at  I  kHz). 

In  a  recent  study  of  attenuation  in  quartz,  sand  in  distilled  water  (which  has  been 
widely  referenced).  Hampton  (Ref.  41)  reported  attenuation  dependent  on  the  square 
root  of  frequency.  The  measurements  do  not  support  this  conclusion.  The  data  (Rcf.4l. 
lig.  ‘))  are  more  in  accord  with  a  first -power  dependency  (as  noted,  also,  by  Mi/ikos  in 
Kef.  50). 

The  values  of  n  which  were  measured,  computed,  and  assumed  above  are  close  to 
that  statistically  computed  by  Attwell  and  kamana  (Ref.  45)  for  frequencies  between  I 
and  I08  Hz,  0.0 1 1 .  Strick  (Ref.  51)  made  a  case  that  the  dependence  of  attenuation 
on  frequency  should  be  close  to.  but  less  than,  I  to  satisfy  causality.  The  writer  believes 
the  best  of  his  experimental  data  are  in  fine  sand  off  Sail  Diego  0/  =  1 .007  ±  0.060). 

When  these  data  are  combined  with  Shumway’s  (Ref.  18)  averaged  data  for  fine  sand. 
ii  =0.W2  ±  0.065. 


Data  listed  in  Tables  I  through  4  and  in  6  are  plotted  in  l  ig.  2  (frequency  versus 
attenuation).  It  can  be  seen  that  most  of  the  data  are  consistent  with  an  approximate 
first-power  dependency  of  attenuation  on  frequency  over  a  wide  frequency  range.  The 
upper  and  lower  bounds  of  the  data  plot  probably  define  the  area  in  which  most  natural 
marine  sediments  will  lie.  With  regard  to  sediment  type,  the  silt-clays  lie  in  a  narrow 
band  at  the  lower  side  of  the  data  plot,  and  very  fine  sands,  silly  sands,  and  sandy  sills 
at  tin.  top.  I  xtrupolution  ot  the  silt-clay  data  to  Ucquc  ikies  below  1  hit/.,  using  a  tilsl- 
power  dependency,  results  in  attenuation  values  in  accord  with  the  data  of  Bennett 
(Ref.  12;  40  to  000  Hz)  and  those  of  Tullos  and  Reid  (Ref.  13:  50  to  400  II/). 

two  re  cc  i  it  studies  of  complex  rigidity  and  c  ml  try  -damping  in  artificial  ciTyvfiarc 
important  implications  in  determining  parameters  for  elastic  and  viscoelastic  models  in 
saturated  sediments.  Cohen  (Ref.  17)  measured  both  p  and  in'  in  complex  rigidity. 

(p  +  in'),  in  artificial,  laboratory  sediments  composed  of  kaolinite  and  bentonite  in  dis¬ 
tilled  water,  with  and  without  a  deflocculating  agent.  Cohen  demonstrated  that  both 
p  and  in'  were  independent  of  frequency  in  the  range  8  6  to  4  7  2  kHz  in  flocculated 
clay;  but  when  a  deflocculating  agent  was  added,  the  flocculated  structure  of  the  clay 
sediment  dispersed,  the  material  lost  all  rigidity,  p.  and  behaved  as  a  Newtonian  fluid  in 
wJWell  theft1  MfMMfti  ditntftttMr  Uf  vum*  i.YM.'Cf’V  U»jU  *"«  taeatfv  dvpxftdottt  up  We- 
quency.  The  addition  ol  35.5  ppt  of  NaCI  caused  reflocculation,  and  complex  rigidity 
was  the  same  as  before.  Kri/.ek  and  Franklin  (Ref.  52).  in  studies  of  shear-wave  energy 
damping  in  flocculated  kaolinite  in  distilled  water,  demonstrated  that  1/(7  was  independent 
of  frequency  in  the  range  0.1  to  30  Hz.  and  that  the  stress-strain  hysteresis  loop  for  a  given 
cycle  was  that  of  a  linear  viscoelastic  medium. 
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The  two  studies  noted  above  (Refs.  17  and  52)  liave  several  important  implications 
for  saturated  clay  sediments,  at  least  in  the  frequency  range  covered  (0.1  !lz  to  43  kHz), 
and  probably  at  much  higher  frequencies:  (1)  saturated,  flocculated  clay  sediments  re¬ 
spond  to  shear-wave  energy  as  linearly  viscoelastic  media;  (3)  the  mu .'pendence  of  energy 
damping  (I  Qs  =  v'  n)  from  frequency  implies  that  linear  attenuation  of  shear  waves 
should  be  proportional  to  the  first  power  of  frequency;  arid  (3)  suspc  is  ions  (without 
flocculated  structures)  do  not  respond  to  wave  energy  as  do  Ho  ecu  la  ted  clay  structures, 
and  almost  all  natural,  high-porosity  silt-clays  have  this  general  type  of  structure. 

Recent  measurements  of  attenuation  in  water-saturated,  natural  sediments  over 
the  frequency  range  3.5  to  1500  k 1 1/  can  be  summarized  as  follows.  Tables  I.  4,  and  b 
list  25  values  of  the  exponent  of  frequency,  n.  between  03>4  and  1 .20.  however,  all  but  two 
of  the  values  fall  between  03)4  and  III.  The  experimental  evidence  indicates  that  the 
dependence  of  attenuation  on  frequency  is  close  toy*,  and  does  not  support  any  theory 
calling  for  a  dependence  of  attenuation  on  f  '1  or  .1  -.  However,  as  in  the  discussion  of 
velocity  dispersion,  the  case  should  not  be  overstated.  As  a  reviewer  pointed  out.  there 
is  no  single  data  set  covering  more  than  two  orders  of  magnitude  in  frequency.  While 
these  data  are  enough  to  show  that  the  dependence  of  attenuation  on  frequency  is  more 
nearly  /  1  than  /  or  /  -,  it  is  not  enough  to  verify  an  exact  dependence. 

REVIEW  OF  ELASTIC  AND  VISCOELASTIC  MODELS 

In  the  field  of  soil  mechanics,  large  static  or  dynamic  stresses  have  to  be  con¬ 
sidered;  and  over  the  full  range  of  stresses,  sediments  may  be  elastic,  viscoelastic,  or 
plastic.  Yong  and  Warkcntin  (Ref.  53.  p.  X0-()4)  present  a  good  discussion  of  the  various 
models  and  elements  within  the  models  which  describe  this  behavior. 

In  the  fields  ol  soil  mechanics  and  foundation  engineering,  the  llookean  model 
and  equations  are  commonly  used  for  derivations  of  dynamic  elastic  constants  and  studies 
of  vibrating  loads  (eg..  Refs.  2X,  31.  32.  and  54-57).  However,  the  dynamic  moduli  from 
most  velocity  data  are  for  very  small  strains,  on  the  order  of  about  I  O'6,  and  corrections 
to  moduli  should  be  made  for  greater  strains  (Ref.  58  presents  a  correction  curve). 

In  the  lields  of  physics  and  geophysics,  studies  of  the  elasticity  of  minerals  and 
rocks  have  demonstrated  that  the  elastic  equations  of  the  llookean  system  adequately 
define  the  velocities  of  compressional  and  shear  waves.  I  hese  equations  are  conveniently 
interrelated  in  a  table  by  Birch  (Ref.  23,  p.  2200).  This  field  has  been  summarized  by 
Birch  (Ref.  50)  and  by  Anderson  and  Lieberman  (Ref.  00).  Papers  of  special  interest 
are  by  Christensen  (Refs.  61, 02).  Brace  (Refs.  03.  ''4).  and  Simmons  and  Brace  (Ref.05). 

The  question  of  water  movement  relative  to  mineral  frame  is  a  critical  key  to 
whether  or  not  the  equations  of  elasticity  can  be  used  in  studies  of  wave  velocities  in 
rock  and  sediments.  If  the  pore  water  does  not  move  significantly  with  respect  to  the 
solids,  then  the  effective  density  of  the  medium  is  the  sum  of  the  mass  of  (lie  water  and 
the  solids  in  a  unit  volume,  there  is  negligible  or  no  velocity  dispersion,  energy  damping 
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is  negligibly  dependent  or  independent  of  frequency,  and  the  equations  of  Hookean 
elasticity  can  be  used  to  study  wave  velocities  and  elastic  constants  within  the  frequency 
range  in  which  these  parameters  are  effective,  unless  attenuation  is  involved  in  the  study. 
This  is  the  "closed  system"  of  Gassmann  (Ref.  66).  The  “closed  system.”  as  a  special 
case  in  studies  of  the  elasticity  or  viscoelasticity  of  saturated  porous  media,  has  been 
noted  in  many  experimental  and  theoretical  studies  (Refs.  I (),  22.  23.  25.  31.  32.  42.  5‘), 
65-81  ). 

Although  the  elastic  equations  of  the  llookean  model  adequately  account  for  wave 
velocities  in  most  earth  materials,  they  do  not  provide  for  wave-energy  losses  in  these 
media.  Fo  account  for  energy  losses,  various  anelastic  (viscoelastic  and  “near-clastic”) 
models  and  equations  have  been  proposed.  Viscoelastic  models  frequently  favored  are  the 
Kelvin-Voigt.  Maxwell,  or  some  other  combination  of  llookean  elastic  springs  and  New¬ 
tonian  dashpots  (see  Ref.  53  for  a  concise  resume),  or  some  variation  of  Biot's  models 
(Ref.  68)  in  which  a  basic  assumption  involves  movement  of  pore  water  of  the  Poise uillc 
type  (at  lower  frequencies). 

In  his  various  papers,  Biot  (eg..  Refs.  68,  61))  discussed  the  full  range  of  systems 
in  which  water  within  pore  spaces  does  or  does  not  move  with  the  solids  upon  imposition 
ol  a  small  stress,  such  as  that  produced  by  a  sound  wave.  In  some  of  these  acoustic 
models,  this  movement  or  flow  of  water  through  the  sediment  mineral  structure  was  con¬ 
sidered  to  be  ol  the  Poiseuille  type.  In  the  last  several  decades,  it  has  been  determined 
that  the  simple  flow  equations  of  the  Poiseuille  type  (derived  from  flow  of  water  through 
tubes)  do  not  hold  lor  real,  in  situ  sediments.  These  equations  have  to  be  consiuerably 
altered,  even  for  clean  sands,  and  are  not  applicable  to  relatively  impermeable  clays  (c.g.. 
Ret.  53).  In  other  words,  models  based  on  Poiseuille-type  flow  of  pore  vater  are  probabh 
not  applicable  to  natural  sediments. 

One  model  which  has  been  especially  studied  iu  connection  with  rocks  and  sedi¬ 
ments  is  the  Kelvin-Voigt  model,  in  which,  as  originally  defined,  eompressional-wave  velocity 
varies  with  frequency;  and  attenuation,  at  frequencies  of  most  interest  in  underwater 
acoustics  and  geophysics,  increases  with  the  square  of  frequency.  White  (Ref.  25. 
p.  I  1 0-1  12)  has  a  thorough  discussion  of  theory  and  experimental  evidence  on  this  sub¬ 
ject  and  concludes  (p.  I  12)  that  neither  velocity  nor  attenuation  shows  this  frequency 
dependence  and  that  the  Voigt  solid  cannot  be  considered  an  adequate  model  of  earth 
materials.  The  evidence  of  this  report  and  earlier  ones  (Refs.  2  and  I'M  are  in  accord 
with  this  conclusion. 


A  VISCOELASTIC  MODEL  FOR  WATER  SATURATED 
SEDIMENTS 

In  the  absence  of  sufficient  experimental  evidence,  it  has  been  possible  to  con¬ 
struct  rather  elegant  theoretical  approaches,  altered— if  necessary— by  constants  to  fit 
available  data.  To  derive  such  theoretical  models,  one  must  start  with  assumptions.  In 
the  case  ol  water-saturated  sediments,  some  ol  the  less  tenable  of  these  assumptions  have 
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Ivon  ill, n  1 1 )  .ill  water-saturated  \oilimonls  arc  .malisons  in  stis|vnsinns  nl  mineral  par- 
Holes  in  fluids:  (2)  all  the  mineral  parlioles  are  spheres:  (3)  Toiseuille  llnw  operates  in 
natural  sediments;  (4)  pore  water  necessarily  moves  relative  to  the  mineral  frame  or  sedi¬ 
ment  structure:  and  (5)  sediments  lack  rigidity,  in  which  case  the  shear  modulus  is  zero, 
and  1‘oisson’s  ratio  is  0.5. 

All  of  the  above  assumptions  are  invalid  in  part  or  in  whole.  I  wo  recent  pa|vrs 
(Refs.  5  and  l‘M  discussed  several  as|vcls  of  sediment  structure  and  elasticity.  Some 
conclusions  are  pertinent,  almost  all  saturated  sediments  have  mineral  particles  which 
are  not  spheres;  near  suspensions  are  unusual;  and  almost  all  sediments  have  sufficient 
rigidity  to  allow  transmission  of  shear  waves.  As  noted  previously.  Vonyi  and  Warkenlin 
(Kef.  53)  indicate  I'oiseuille  flow  4 1 hroiiyrh  small  tubes)  does  not  hold  lor  natural  sediments 

(liven  macroscopic  isotropy,  small,  sinusoidal  stresses,  wavelengths  much  (treater 
than  (train  si/e.  and  frequencies  from  a  few  hertz  to  at  least  several  hundred  kilohertz 
(and  probably  into  the  megahertz  range  for  most  natural  sediments),  some  parameters 
in  addition  to  those*  in  the  preceding  paragraph  are  as  noted  m  previous  sections:  attenua¬ 
tion  in  decibels  per  linear  measure  is  approximately  dependent  on  the  first  power  of  Ire- 
tpicncy :  and  velocity  dispersion,  if  present,  is  small.  Some  relative  movement  of  pore- 
water  and  mineial  frame  cannot  be  excluded  on  the  basis  of  present  evidence,  although 
the  above  parameters  indicate  that,  if  present,  it  should  be  small. 

I  he  model  proposed  below  is  intended  as  a  tentative,  working  model.  It  should 
be  emphasized  that  other  models  are  not  excluded  if  they  are  within  the  stated  parame¬ 
ters.  Ihe  whole  subject  merits  much  more  experimental  and  theoretical  study. 

A  model  and  concomitant  equations  within  Ihe  parameters  noted  above  is  a  case 
of  linear  viscoelasticity.  The  basic  equations  of  linear  viscoelasticity  have  been  summarized 
in  an  excellent  treatise*  by  l  erry  (Ref.  82).  l  or  the  model  recommended  in  this  paper, 
the  basic  equations  (Ref.  83)  have  been  discussed  in  different  form,  including  neglect  of 
negligible  factors,  in  Ref*  2.  25.  52.  82  (pp.  ‘>3.  4).  84  and  85. 

In  the  above  model.  Ihe  fame  elastic  moduli  /a  and  X  are  replaced  by  complex 
moduli,  (p  +  ip')  and  (X  +  /X').  in  which  p,  X,  and  density  govern  wave  velocity,  and  Ihe 
imaginary  moduli,  ip  and  /X'. govern  energy  damping.  The  following  (Ref.  82.  pp.  1 1*13) 
illustrate  the  stress-strain  relations  in  this  model.  Tor  a  sinusoidal  wave,  if  the  viscoelas¬ 
tic  behavior  is  linear,  the  strain  will  be  out  of  phase  with  stress.  I  he  stress  can  be  vec- 
tonally  decomposed  into  two  components:  one  in  phase  with  strain  and  one  ‘>0  deg  out 
of  phase.  Tor  a  shear  wave,  the  complex  stress/strain  ratio  is  p*  -  p  +  ip'.  The  phase 
angle.  <>.  which  expresses  energy  damping  is,  in  this  case*:  Ian  0  =  p'/p. 

The  basic  derivations  of  the  above  model  are  in  Refs.  25  and  82  and  will  not  be 
repeated  here.  Without  assumptions  as  to  negligible  factors,  the  equations  of  the  model 
in  the  form  used  by  Rucker  ( in  Ref.  2.  p.  404b)  or  by  Terry  ( Ref.  82.  p.  04.  410).  reduce 
to  the  following  for  both  compressional  and  shear  waves  (with  some  changes  in  notation). 
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I  (J  is  I  lie  specific  attenuation  factor  or  specif:.  dissipation  function 
ii  is  I lu*  attenuation  coefficient 
I"  is  wave  velocity 

I  is  Ireiiiieney  (circular  frequency.  o :  2nb 

Subscripts  (/*or  o  can  Iv  inserted  into  I  <|  (2>  when  referrint'  (o  comprcssional 
or  shear  waves, 

When  energy  dampum  is  small  (/  <■  .  X’  •  * .  X  and  p'  •  '<  //:  Kel.  25.  p  ‘>5.  Ret.  s. 
p.  125  r  «  I.  where  r  */ 1 '  2irl\.  the  term  in  the  denominator  of  I  «|.  ( 2 ».  ./*  I  •’ 
is  ncy'liyihlc  and  can  be  dropped.  Mils  leaves  the  more  familiar  expression  (re  .  Rets  2! 
42.  44.  and  45 1. 
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Where  (in  addition  to  those  symbols  alreatlv  defined  I 


A  is  the  logarithmic  decrement  (lop  of  the  r.ilio  of  two  successive  amplitudes  m 
.m  exponentially  decaying  siiuiviiil.il  w.oe) 

tan  0  is  the  loss  angle 

A/'.  /.  is  tract i« in  of  strain  energy  lost  per  siress  cycle 

(t  is  attenuation  in  tlecilvls  per  linear  measure  <«  g  .  dll  cm) 

I  ipi.il ions  involving  eompressional-  and  shear-wave  velocities  in  Kel.  2.  or  in 
Kef.  7<>  are  tin  I  erry's  notation): 

(X  +  2p)  -  pl,;t  I  -  r*  )/t  I  +  r-  )- 

P  I  -/-Ml  +  /  - I- 

wlierc 


t‘>) 

(10) 


r-  ,i  I  2iri 
X  r  I  . line’s  constant 
p  rigidity 
fj  -  density 

In  I  (|\.  ('))  ami  ( 10).  the  term  ( I  -  /;)  ( I  +  r- )-  indicates  the  degree  of  velocity 
dispersion  lor  linear  viscoelastic  media.  When  damping  is  small  (defined  above),  this  term 
is  negligible  and  can  be  dropped,  as  implied  in  Kef.  X2.  p.  04  This  leaves  the  more 
familiar  llookean  equations. 


(X  f  2p)  pi,; 

ill) 

p  -  p»;J 

(12) 

I  his  means  that  if  the  factor!  I  -  r' )  <  I  +  r- )-’  in  I  qs.  (0)  and  (10)  is  considered 
negligible  and  dropped,  wave  velocity.  I !(].  and  the  log  decrement  are  independent  of 
frequency.  and  linear  attenuation  is  proportional  to  the  first  power  of  frequency. 

Computations  performed  with  the  data  of  this  report  and  from  the  literature  indi¬ 
cate  that  most  water-saturated  rocks  and  sediments  qualify  under  the  above  definitions  as 
media  with  “small  damping.”  Therefore.  I  qs.  t.h.  (4).  (II).  and  ( 12)  should  apply  to  both 
water-saturated  sediments  and  rocks.  However,  those  investigators  who  wish  to  include 
velocity  dispersion  and  I  <J  or  a  log  decrement  dependent  on  frequency  can  consider 
I'.qs.  (2).  O).  and  ( 10). 
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In  Ihi'  sediments  discussed  in  this  report.  t Ik*  quality  factor.  nr  specific  dissipation 
function.  (J.  I  (J.  and  (lie  logarithmic  decrement.  2k.  arc  approximately  independent  of 
frequency.  Hus  follows  I'ront  t Ik*  dependence  of  attenuation  (approximately  I  on  the 
first  power  of  frequency.  Table  7  lists  these  properties  for  the  hi  \ilii  measurements  in 
Tables  I  and  2.  A  conclusion  from  llu  data  nl  this  report  is  that  the  approximate  fre¬ 
quency  independence  of  I  IQ.  (J.  and  ^  can  be  extended  from  rock  (cj a .  Ref.  4.*  I  through 
most  natural,  water-saturated  sediments. 
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Jlld  Ilk'  l»)!JllllllMI.  IkSH'IIU'lll.  Jt  .  I  loin  In  Sun  Measiiieiiuiils<i| 
<nni|>tessitinal\V.ise  Veins  it).  I„.  and  Attenuation.  in  Sediments 
Off  San 


Sedinienl  1  s  pe 

tv 

•  tv 

s>' 

Nules 

Si  lid 

Coarse 

32 

<1.1131 

0.011 

1 

21 

0.034 

0.107 

2 

Medium 

31 

0.032 

0.101 

1 

35 

0,028 

0.081 

s 

1  me 

31 

0.032 

0.100 

1 

31 

0.032 

(Mill 

s 

44 

0.023 

0.071 

3 

Very  fine 

32 

0.031 

0.013 

1 

21 

0.034 

0.107 

> 

Sandy  sill 

23 

0.044 

0.138 

s 

54 

0.018 

0.058 

3 

Sand-silt-elay 

31 

0.033 

0.102 

s 

Oayey  sill 

104 

0.010 

0.030 

s 

III 

0.001 

0.028 

1 

114 

0.001 

0.028 

1 

IIX 

0.001 

0.027 

3 

21.3 

0.004 

0.012 

1 

437 

0.002 

0.007 

1 

Sand-sili-elay 

368 

0.003 

0.001 

1 

*  Computed  from  e(|u.itinns 


I IQp  =  “p I ‘pM  attenuation,  «/(  H  f> Kf*  ,1^ 


Nolos-  |.  Table  I  (fan  A I;  frequency :  14  kll/. 
2.  I'ublc*  2;  frequency:  Kill  kll/. 

.3.  table  I  (fart  H);  frequency:  25  kll/. 
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PROPIRTIIS 


Introduction 

The  relationships  between  attenuation  anil  other  physical  properties  in  saturated 
sediments  are  of  considerable  importance  in  determinini!  the  causes  of  attenuation  and 
in  selecting  appropriate  anelaslic  models.  In  this  section,  various  relationships  will  be 
briefly  noted  and  illustrated  prior  to  discussions  of  the  causes  of  attenuation. 

As  discussed,  attenuation  in  decibels  per  linear  measure  is  approximately  dependent 
on  the  first  power  of  freipiency :  that  is.  in  the  equation  a  -  A/".  n  is  close  to  I.  If  n  is 
taken  as  I.  the  only  variable  in  the  equal  ions  lor  various  sediments  is  the  constant  A.  This 
constant  is  particular!)  useful  in  relating  attenuation  toother  sediment  properties,  such 
as  grain  si/e  and  porosity.  The  relations  between  A  and  other  physical  properties  give  an 
insight  into  the  causes  of  attenuation,  which  allows  prediction  of  attenuation  (as  disit.  .xcd 
in  the  last  section)  after  deriving  a  value  of  A  from  its  relationships  with  mean  grain  si/e 
and  porosity. 

With  an  assumption  that  linear  attenuation  is  dependent  on  the  first  power  of  fre¬ 
quency.  values  of  A  can  be  easily  computed  by  dividing  attenuation  by  frequency.  This  was 
done  for  all  data  in  Tables  I  to  4  except  as  follows:  ( I  I  the  values  of  A  from  Table  I  (I’art  A) 
were  determined  from  the  14-kll/  measurements  (considered  most  reliable):  and  (2)  for 
three  measurements  in  Table  4a  which  did  not  show  a  first-power  dependency.  A  was  de¬ 
termined  from  an  attenuation  computed  at  a  frequency  near  the  mid-range  of  the  mea¬ 
surements.  These  values  of  A  were  then  plotted  versus  mean  grain  si/e  and  porosity  in 
figs.  3.  4.  and  5. 

Relationships  between  A  (or  attenuation)  and  mean  grain  si/e.  porosity,  and  dy¬ 
namic  rigidity  are  illustrated  and  discussed  in  the  following  sections.  The  data  are  for 
natural  sediments:  both  in  situ  and  laboratory  measurements  are  included.  Averaged 
values  from  Tables  I  through  3  were  used  as  primary  data  to  establish  regression  equa¬ 
tions  (in  the  figure  captions),  but  the  individual  measurements  which  were  averaged  are 
shown  to  better  illustrate  the  trends  and  scatter  of  the  data.  These  regression  equations 
arc  included  for  use  in  predicting  attenuation  when  grain  size  and  porosity  are  known. 
These  equations  are  strictly  empirical  and  are  recommended  only  within  the  limiting 
values  indicated.  The  values  of  A  so  obtained  are  approximations,  but  it  is  predicted  that 
most  future  measurements  of  attenuation  will  result  in  A  values  which  fall  within  the  indi¬ 
cated  “envelopes.” 

Attenuation  and  Grain  Si/e 

When  wavelengths  approach  grain  si/e.  Rayleigh  scattering  takes  place,  and  attenua¬ 
tion  is  proportional  to  the  fourth  power  of  frequency.  This  effect  has  been  discussed  by 
lCik*|Kiff  nod  Potter  (Ref  in  tfvc  cam*  gt«H)(e.  RayleigU  rovfWwtg  iw»*Ht*  U*  K*  a 
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Mean  Grain  Si/e 

I  igurc  .V  Mean  grain  size,  M,,  in  e>  uml,  ami  microns,  u.  versus  k  (in  u  k/' )  in  natural, 
saturated  sediments  ( tables  1  -4a );  see  test.  Solid  symbols  are  averaged  values  lor  data 
oil  San  Diego:  iliamonds- 1 9<»6- 1 4 7(1 ;  lrumglcs-|956;  inverted  lri.mglcs-1962:  eireles- 
Shutiivvav  (Kef,  ISl,  Table  5;  open  svmbols-vlata  in  averages;  X  literature  value  from 
Table  4a.  Area  between  dashed  lines:  predicted  area  within  which  most  data  should 
fall.  Regression  equations  for  solid  lines  (for  data  off  San  Diego  and  selected  literature 
values-reconunended  for  use  in  similar  sediments):  Coarse,  medium,  and  fine  sand,  in 
part  (0  to  2.6  <•):  k  0.4556  +  0.0245(4/.).  I  me  sand  on  part ).  very  fine  sand,  and  mixed 
sizes  ( 2.6  to  4.5 0):  k  11.1978  +  0.12450/.).  Mixed  sizes  (4.5  to  6.0 <,v):  k  ■  8.0599  - 
2.52280/.)  +  0.200l)8  l.l/-)2.  Sill-clays  (6.1?  to  9.5  o).  k  =  0.9451  •  0.20410/.)  +  0.01 170/.I2. 


(actor  in  Bushy  ami  Richardson's  (Ret.  21  )  measurements  ol  attenuation  in  the  megahertz, 
range  in  sands.  I ‘or  most  sediments,  Rayleigh  scattering  is  not  a  factor  to  at  least  several 
hundred  kilohertz,  and  probably  into  the  megahertz  range. 

f  igure  d  illustrates  the  relationships  between  the  constant  A  (in  a  hi 1 )  and  mean 
grain  size  for  most  of  the  sediments  listed  in  Tables  I  to  4.  Mean  grain  si/e  is  plotted  in 
logarithmic  phi  units  (0  =  -log,  of  grain  diameter  in  millimeters;  see  Ref.  N7  for  discus¬ 
sion).  In  studying  f  ig.  4,  one  can  easily  translate  k  into  linear  attenuation  in  decibels  per 
meter,  because  at  any  given  frequency  the  only  variable  is  A 
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porosity  mixed  sizes  <46.7  to  52  percent  I  A  -  0.0490300  -  1.7688.  Mixed  sizes  <52 
to  65  percent):  A  =  3.3232  -  0.0489<n).  Silt-clays  <65  to  90  percent):  A  =  0.7602  - 
0.01487OD  +  0.0000780D2. 


Figure  3  illustrates  the  distinct  difference  between  mean  grain  size.  4/,.  and  the 
constant  A  (or  attenuation)  in  sands  (,!/.  =  0  to  4  0,  or  I  to  0.0625  mm)  and  in  the  silt- 
clays  (0  >  4),  In  sands  d  igs.  3  and  4).  A  increases  gradually  with  0  (or  decreasing  grain 
size)  from  coarse  into  fine  sands  (to  about  2.6  0),  and  then  A  increases  rapidly  into  the 
finer  sand  sizes.  Other  investigators  have  measured  increased  attenuation  in  sands  with 
decreasing  grain  size  (e  g..  Kel's.  8,  18.  and  41.  I;ig.  ')).  The  maximum  values  of  A  are  in 
very  fine  sand  and  in  mixtures  of  sand,  coarse  silt,  and  clay  (<\g..  silty  sand,  sandy  silt) 
in  the  grain-size  range  of  3.5  to  4.5  0  (0.0‘>  to  0.04  mm).  Fquivalcnt  lower  values  of  A 
are  in  coarser  grained  sand,  in  silt,  and  in  finer  sized  sandy  silt.  Attenuation  decreases 
with  decreasing  grain  size  (increasing  0)  from  about  4.5  0  to  about  6  0. and  then  gradually 
declines  with  decreasing  grain  size  into  fine  silts  and  clays. 

A  semi-log  plot  (using  phi  units:  Fig.  3)  tends  to  clarify  relationships  between 
grain  size  and  other  properties  in  the  finer  sizes  (silt-clays),  but  obscures  relationships  in 
sands.  Therefore,  mean  grain  size  in  millimeters  is  plotted  versus  A  in  Fig.  4.  Figure  4 
illustrates  better  than  Fig.  3  the  gradual  increase  of  A  from  coarse  and  medium  sand  sizes 
into  fine  sands,  and  the  marked  increase  in  A  from  about  0.17  mm  (2.6  0)  into  the  finer 
sand  sizes. 

Attenuation  and  Porosity 

Porosity  varies  with  A  (Fig.  5)  in  the  same  way  that  mean  grain  size  varies  with  A 
(Figs.  3  and  4).  Fquivalcnt  values  of  A  are  apt  to  be  found  in  the  coarser  sizes  of  sand 
and  in  higher  porosity  silts  and  sandy  silts.  The  highest  values  of  A  are  in  very  fine  sands 
silty  sands  and  sandy  silts  in  the  porosity  range  ol  50  to  54  percent.  In  silt-clays.  A  de¬ 
creases  with  increasing  porosity.  Sluimway  (Kef.  18.  Fig.  6)  indicates  a  similar  variation 
tv t ww«  nUi'nuaiioo  mul  pmntiiy  ai  CtvquvJJcie*  Cviwwn  v  -uul  Ju  UU 
between  grain  size  and  porosity,  as  they  affect  attenuation,  will  be  discussed  in  a  section 
below. 

Attenuation  and  Dynamic  Rigidity 

Figures  6  and  7  (reproduced  from  Kef.  10)  illustrate  the  probable  variations  of 
computed  values  of  dynamic  rigidity,  /a,  with  grain  size  and  porosity.  C’onr  arison  of 
these  figures  with  Figs.  3,  4.  and  5  indicates  that  A  and  values  of  rigidity  re*  v'n<l  in  the 
same  way  to  variations  in  grain  size  and  porosity. 


PROPOSED  PREDICTION  METHOD 

CAUSES  OF  ATTENUATION 
Introduction 

In  studies  of  the  attenuation  of  compressional  and  shear  waves  in  saturated  porous 
media,  there  are  usually  two  (or  a  combination  of  two)  common  viewpoints:  (I )  a  visco¬ 
elastic  model  is  used  in  which  pore  water  moves  relative  to  mineral  grains,  and  pore-water 


figure  (>.  Moan  grain  si/c  versus  computed  values  of  the  shear  modulus  ulvuaimc 
rigid  it  v);  from  llamiltou  (Kef.  19).  Square^  are  mean  values;  crosses  are  .'  times 
the  standard  error  of  the  mean.  T-samples  from  the  continental  terrace  (shelf 
and  slope):  I’-alnssal  plain:  ll-ahvssal  hill  lower  left:  St.  IV ter \  sand  (circles ) 
and  Ottawa  vmd  (triangles).  See  Ref.  19  for  discussion. 

viscosity  and  media  permeability  are  dominant  factors  in  viscous  sound  absorption:  or 
(2)  a  linear  viscoelastic,  or  “nearly  elastic,"  model  is  favored  in  which  attenuation  is 
mostly,  or  entirely,  due  to  internal  friction  energy  is  lost  in  intercrystalline  or  inter- 
grain  movements). 

The  case  has  been  made  that  in  crustal  rucks,  both  dry  and  saturated,  internal 
friction  is  by  far  the  most  probable  dominant  process  in  wave-energy  damping  at  seismic 
frequencies  to  at  least  several  hundred  kilohertz,  and  that  viscous  losses  are  probably 
negligible  (c.g.,  Refs.  25  and  88).  The  writer  believes,  as  discussed  below,  that  internal 
friction  is  by  far  the  dominant  dissipative  process  in  water-saturated  sediments. 

It  is  apparent  from  the  present  study  that  sands  with  grain-to-grain  contacts  and 
no  cohesion  (physicochemical  net  attraction)  should  be  studied  separately  from  silt-clays 
with  cohesion  and  mineral  particles  probably  separated  by  layers  of  adsorbed  water. 
Consequently,  the  following  discussion  will  deal  separately  with  these  two  structural 
types.  The  structures  of  these  and  other  sediment  types  were  reviewed  in  Refs.  5  and  l‘>. 
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Porosity,  percent 

I  iguro  7.  Sediment  porosity  versus  dynamie  rigidity;  symbols  and  remarks  as 
in  fig.  6 


Attenuation  in  Sand 

Study  of  attenuation  in  sands  is  of  particular  interest  in  determining  the  causes  of 
attenuation  because  these  sediments  are  the  most  permeable,  and  viscous  losses,  if  present 
should  be  higher  in  sand  than  in  other  sediment  types.  On  the  other  hand,  if  frictional 
energy  losses  are  dominant,  the  intergrain  reactions  to  stress  which  cause  attenuation 
should  he  related  to  the  same  complex  factors  which  affect  dynamic  rigidity  and  static 
shear  strength.  These  factors  will  be  briefly  reviewed.  Reference  l‘>.  a  more  extensive 
review,  contains  numerous  references  to  items  outlined  below. 

Sands  have  grain-to-grain  contacts  between  mineral  particles,  and  resistance  to 
shear  stress  is  related  to  sliding  and  rolling  friction  between  grains,  to  the  number  of 
intergrain  contacts,  and  interlocking  between  grains.  See  Kef.  53  for  a  good  resume 

I  lie  number  of  interpart  iele  contacts  m  sands  depends  on  grain  si/e  and  density  of 
packing  1/  r  .  loose,  dense,  or— most  frequently— an  intermediate  packing!.  At  any  giyen 
particle  si/e.  porosity  is  a  measure  of  packing  Interlocking  of  grains  increases  with 
density  of  packing  and  angularity  of  grains. 


Laboratory  experiments  with  clean  quartz  sand  grains  (St.  Peter's  or  Ottawa  sand) 
have  indicated  the  following  (Kef.  28):  (I)  at  the  same  grain  si/e.  dynamic  rigidity  in¬ 
creases  with  decreasing  porosity  in  sands  with  round  grains  because  of  more  interparticle 
contacts  with  denser  packing;  (2)  at  the  same  grain  si/e  and  porosity,  sands  with  angular 
grains  have  higher  dynamic  rigidities  and  shear-wave  velocities  than  sands  with  round 
grains  because  of  interlocking  between  grains:  and  (3)  at  the  same  grain  si/e.  increased 
angularity  of  grains  (causing  greater  rigidity)  can  be  more  important  than  increases  in 
porosity  (causing  lesser  rigidity). 

The  grains  of  natural  marine  sands  are  much  more  angular  than  those  in  St.  Peter’s 
or  Ottawa  sand  at  any  given  grain  si/e.  In  natural  sands,  angularity  increases  and  sphericity 
decreases  with  decreasing  grain  size  (Ref.  8{)).  Interlocking  of  grains  should  be  greater 
when  coarse  silt  particles  are  present. 

The  discussion  above  and  the  conclusions  presented  in  Ref.  I()  imply  and  predict 
for  natural  marine  sediments; 

(1 )  Coarse  sand  has  fewer  intergrain  contacts  and  the  grains  are  rounder  than  in 
the  finer  sizes:  thus,  dynamic  rigidity  should  be  a  minimum  in  coarse  sand 
because  of  lesser  intergrain  friction. 

(2)  Dynamic  rigidity  should  be  less  in  sands  composed  of  the  highly  rounded 
grains  of  St.  Peter's  or  Ottawa  sand  than  in  natural  sands  of  the  same 
density,  porosity,  and  grain  si/e:  thus  tests  with  St.  Peter's  or  Ottawa  sand 
cannot  always  be  directly  related  to  natural  sands. 

(3)  Hardin  and  Richart  (Ref.  28)  observed  that  an  increase  in  porosity  will  cause 
a  decrease  in  rigidity.  This  is  only  true  in  sands  of  the  same  grain  size  and 
angularity. 

(4)  In  natural  sands,  porosity  increases  with  finer  si/es  which  have  more  numerous 
intergrain  contacts  and  more  angular  grains,  which  may  cause  rigidity  to  in¬ 
crease  with  increasing  porosity. 

(5)  At  some  relatively  fine  sand-silt  si/e.  intergrain  friction  and  interlocking  will 
reach  a  maximum  and  rigidity  will  be  at  a  high  point  for  natural,  uncemented 
sediments. 

((>(  When  sand  grains  are  no  longer  in  contact  because  of  increased  amounts  of 
fine  silt  and  clay,  intergrain  friction  and  interlocking  are  no  longer  relatively 
effective  and  rigidity  is  mostly  due  to  cohesion  between  finer  particles. 

(7)  Cohesion  between  fine  particles  varies  with  the  numerous  factors  discussed 
below:  but.  in  general,  cohesion  and  rigidity  decrease  with  increasing  porosity 
and  decreasing  grain  si/e  in  higher  porosity  silt-clays. 


Figures  b  ami  7  illustrate  Ilk  phenomena  outlined  above.  Of  the  sands,  coarse 
sand  lias  the  least  rigidilv.  and  natural  sands  have  higher  rigidities  than  round -grained 
St  Peter's  sand.  Rigidilv  increase's  slurplv  with  decreasing  train  si/c  and  increasing 
porosity  to  maximum  values  in  very  line  vuuls.  silty  voids,  and  coarse-  silts.  I  liese  ina\i- 
muni  values  of  dynamic  rigidity  occur  between  mean  train  si/es  of  4.5  to  4.5  <?  lO.O*)  to 
<'.04  mill)  and  porosities  between  >5  and  <»()  percent.  In  a  density  veTsiis  dvnaniic 
rigidity  plot  (not  shown*,  ilia \111uun  rigidities  occur  between  densities  of  1.7  and  I  .X 
g  cm’.  As  trains  become  finer,  rigidity  decreases  with  increasing  porosity  and  decreasing 
density  (l  it.  7)  and  decreasing  grain  si/e  (I  ig.  h).  Scatter  diagrams  of  percent  sand  si/e 
by  weight  versus  rigidity  (not  shown)  indicate  nmiiinim  rigidity  normally  occurs  when  a 
sediment  is  composed  of  about  (>0  to  r»5  percent  void  In  weight. 

Dynamic  rigidity  is  a  measure  of  resistance  (friction  and  interlocking  between 
grains)  to  shearing  forces,  which  tend  to  move  grains.  If  attenuation  is  due  to  energy  lost 
by  friction  between  grains,  then  rigidity  and  attenuation  should  vary  because  of  the  same 
factors.  However,  as  with  rigidity,  there  is  complex  interaction  between  grain  si/es.  shapes, 
density  of  packing,  or  porosity. 

Mean  grain  size  is  an  important  factor  in  both  rigidity  and  attenuation.  In  coarse 
sands,  the  grains  are  larger  and  more  rounded,  interparticle  contacts  fewer,  and  surface 
areas  smaller,  than  in  the  finer  sands.  Consequently,  rigidity  is  relatively  low  in  coarser 
grained  sands;  and  when  grains  are  moved  the  fewer  interparticle  contacts  produce  only 
a  small  attenuation  that  can  be  attributed  to  friction  (bigs.  .7  and  4).  As  grain  size  de¬ 
creases,  the  grains  become  more  numerous  in  a  unit  volume,  they  are  more  angular,  and 
consequently,  have  larger  surface  areas  and  more  interparticle  contacts.  Therefore,  rigidity 
increases.  When  grains  are  relatively  moved,  however,  the  more  numerous  interparticle 
contacts  and  greater  surface  areas  result  in  greater  energy  losses.  This  has  been  shown  in 
the  laboratory  (rg.,  Kef.  41.  Fig.  <))  and  in  situ  (this  report  and  Ref.  X). 

flic  change  in  rate  of  increase  of  A  from  coarser  sands  at  about  2.U  <?.  or  0.1  7  mm. 
into  the  finer  sand  si/es  (Figs.  4.  4)  requires  further  discussion.  The  change  is  apparently 
related  to  the  relationships  between  grain  si/c.  porosity,  and  surface  area,  or  the  number 
of  interparticle  contacts. 

A  semi-log  plot  of  mean  grain  si/e  (in  phi  units)  versus  porosity  (Kef.  7.  l  ig.  2) 
indicates  an  almost  linear  relationship.  This  type  of  plot  obscures  the  real  relationship. 
Figure  X.  an  arithmetic  plot  of  grain  si/c  (in  millimeters)  versus  porosity.  illustrates  the 
very  gradual  increase  of  porosity  with  decreasing  grain  size  in  the  coarse  and  medium 
sand  si/es.  and  the  marked  increase  in  porosity  which  occurs  as  grain  si/es  decrease  from 
fine  sands  into  the  silt-day  si/.cs  (porosities  greater  than  about  b()  percent).  There  is 
much  scatter  in  the  relationships  between  grain  si/e  and  porosity,  because  at  any  given 
grain  si/e.  porosity  varies  with  packing  (/.<  .,  dense,  loose). 

II  attenuation  is  related  to  energy  lost  in  intergrain  friction,  a  critical  factor  in¬ 
volves  the  number  of  intergrain  contacts,  or  surface  areas  involved,  per  unit  volume.  If 


Porosity,  percent 


I  ii*iii s'  x.  Si’tliiiH'til  pi»rosit\  u'imis  tiHMit  L'r.nn  si/c  in  imn.  SvinhoK  ilumoiuls- 
frmn  l.ihK's  I  .pul  iin  K's- l.il'lo  horn  Rot .  I  K;  M|ii.in-s-fioiii  Kit 

grain  si/es  decreased  without  a  change  in  porosity.  the  more  numerous  contacts  should 
lead  to  ,m  increase  in  attenuation.  Il  porosity  increaseil  without  .1  chance  in  .main  si/e. 
the  fewer  inlergrain  contacts  should  lead  to  decreased  attenuation.  In  coarse  and  medium 
sands  there  is  a  relatively  small  increase  in  porosity  with  decrease  in  main  si/e  (Tig.  K). 

As  grain  si/es  decrease  from  0.500  to  0.1  SO  mm  (I  to  2.5 v>>.  there  is  an  increase  in 
porosity  of  about  5  to  (>  percent,  and.  accord iii{!  to  Slmmway  and  Igclman  (Kef.  ‘>0 >. 
grain  surface  areas  increase  by  about  240  percent,  from  0.1  SI)  to  0.0(>5  mm  (2.5  to4</0. 
porosity  increases  about  (>  to  10  percent,  but  grain  surface  area  increase’s  In  about  520 
percent.  In  summary,  as  grain  si/es  decrease  from  0.500  to  0.0(0  mm  1 1  to  4  0 1 1  here  is 
an  increase  of  about  15  percent  porosity,  but  grain  surface  areas  increase  on  the  order  of 
bOO  to  700  percent.  If  attenuation  is  due  to  energy  lost  between  mineral  particles  in 
contact,  then  the  increase  in  surface  area  (more  mineral  particles  in  contact  in  a  unit 
volume)  should  cause  increased  attenuation  (A  higher)  which  is  not  offset  by  greater 
porosity.  In  addition,  the  rate  of  increase  in  attenuation  should  be  less  from  0.500  to 
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0.180  mm  ( coarse-mod iuin  sand  into  lino  sand)  than  (ho  rate  from  about  0.1X0  mm  into 
fill  r  sizes  (the  increased  rate  actually  occurs  at  about  O.l-7  nun:  Figs.  3.  4). 

Figure  5  i  porosity  versus  k  »  illustrates  these  same  relationships:  k  increases  gradu- 
ly  as  porosity  unreases  from  5.X  to  about  47  percent,  and  then  increases  sharply  from 
is  porosity  range  to  porosities  of  about  52  percent.  These  increases  in  attenuation  with 
.ncreasing  porosity  in  sands  are  related  through  grain  size  and  the  more  numerous  par¬ 
ticles  and  greater  grain  surface  areas,  as  discussed  above. 

Attenuation  in  Fine  Silts  and  Clays 

In  fine  silts  and  clays,  shear  strength,  dynamic  rigidity,  and  attenuation  are  ap¬ 
parently  related  to  cohesion  between  fine  particles.  Cohesion  is  the  resistance  to  shear 
stress  which  can  be  mobilized  between  adjacent  fine  particles  which  stick,  or  cohere,  to 
each  other.  Cohesion  is  considered  to  be  an  inherent  property  of  fine-grained,  clayey 
sediments  which  is  independent  of  stress;  it  is  caused  by  physicochemical  forces  of  an 
interparticle,  intcrmolecular.  and  intergranular  nature.  This  subject  is  reviewed  in  Ref.  DC 
which  contains  numerous  references  to  the  following  items.  Clay  particles  may  not  be  in 
contact;  they  are  apparently  surrounded  by  layers  of  adsorbed  water  through  which  they 
interact  with  other  particles.  The  amount  of  pore  water,  the  distance  between  particles, 
and  the  number  of  interparticle  contacts  are  important  influences  on  cohesion.  At  points 
of  near  contact  between  clay  particles,  there  is  often  bonding  akin  to  cementation,  espe¬ 
cially  in  the  presence  of  iron  oxides,  calcium,  silica,  and  other  minerals  in  solution  in 
interstitial  waters.  Where  sediments  have  been  exposed  to  overburden  pressures,  there  is 
apt  to  be  pressure-point  solution  and  redeposition.  The  structure  of  the  mass  of  clay 
particles  is  important;  for  example,  it  has  been  demonstrated  that  the  flocculated,  or 
“eardhouse”  structure  (Ref.  5.  Fig.  Id)  is  the  strongest.  These  structures  are  largely  de¬ 
termined  In  interparticle  forces  and  the  number  of  interparticle  contacts.  Differing  day 
minerals  affect  cohesion  because  of  particle  size  and  differing  interparticle  forces;  for 
example.  Na-monlmorillonite  has  stronger  cohesive  bonds  than  kaolinite.  Shear  stress  in 
clayey  sediments  occurs  bet  we  n  particles  and  not  through  them,  near-contact  points  will 
deform  elastically  or  plastically  (depending  on  stress)  by  an  amount  to  sustain  the  effec¬ 
tive  stress.  Homogeneous  clays  are  practically  impervious;  for  example,  the  coefficient 
of  permeability  (centimeters  per  second)  in  clean  sand  is  of  the  order  of  10  -\  whereas 
this  coefficient  in  homogeneous  clay  is  of  the  order  of  10  K  (Ref.  l>|  ). 

The  highest  attenuation  in  marine  sediments  is  in  very  fine  sand  and  coarse  silt 
and  in  mixtures  of  sand  ami  silt.  As  grain  sizes  become  smaller,  the  sand  and  coarser  silt 
grains  become  fewer  and  separated,  and  relative  movement  between  grains  is  controlled 
by  complex  interparticle  forces  (cohesion).  Cohesion,  dynamic  rigidity .  and  attenuation 
become  less  as  grain  sizes  decrease  and  porosities  increase  to  the  level  of  hiuh-porosity 
silt-days.  These  effects  are  apparently  due  to  the  same  causes:  the  strength  of  interpar- 
ticle.  net  attractive  forces  and  the  number  of  interparticle  contacts.  These  forces  become 
weaker  and  contacts  fewer  with  higher  porosities,  given  the  same  mineral  and  structural 
type. 
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In  (lie  curves  relating  A  to  porosity,  the  slurp  inlleetion  wlneh  occurs  jt  about 
(»5  percent  tl'ig.  5l  and  at  a  mean  grain  si/e  o!  about  (<  0  40  01  «*  nun:  lags.  4.  4)  indi¬ 
cates  the  probable  range  of  porosity  and  mean  grain  si/e  where  the  separation  of  larger 
grains  (sand  and  silt  I  occurs.  McCann  and  McCann  (Ref.  X>  also  selected  (>  0  as  the  gram 
si/e  at  which  interparticle  forces  become  dominant. 

It  is  important  in  studies  of  attenuation  in  clay-water  sy stems  to  understand  the 
different  mineral  structures.  Clay,  in  the  presence  of  an  electrolyte  (e  g  .  seawater!. forms 
flocculated  structures  which  have  a  finite  rigidity  and  transmit  shear  waves  (see  Ref.  Id). 
When  this  flocculated  structure  is  dispersed  and  the  physicochemical  c<*uds  are  broken, 
the  mixture  is  an  actual  or  near  suspension.  Attenuation  in  dispersed  or  flocculated  clay- 
water  systems  is  distinctly  different.  Cohen  (Ref.  17)  demonstrated  that  both  p  and  ip 
were  essentially  independent  of  frequency  in  flocculated  kaolinitc  in  distilled  water;  but 
when  a  defloceulating  agent  was  added,  the  mixture  lost  rigidity,  g.  and  behaved  as  a 
Newtonian  fluid.  This  means  that  I  IQ  for  shear  stresses  (p'  p)was  independent  of  fre¬ 
quency.  and  attenuation  dependent  on  the  first  power  of  frequency  in  the  flocculated 
structure;  but  in  the  dispersed  mixture,  attenuation  was  dependent  on  ./-.  When  Cohen 
added  NaCI.  the  mixture  reflocculated  and  complex  rigidity  was  the  same  as  before. 

I)e  Cl  raft -Johnson  (Ref. ‘>2)  experimentally  showed  that  dispersed  clay  structures  have 
higher  log  decrements  than  flocculated  structures  at  the  same  water  content.  McCann 
(Ref.  (>3)  explained  this  effect  theoretically.  T.xcept  as  an  approach  to  understanding 
theory,  experimental  data  on  attenuation  in  dispersed  day-water  systems  are  not  appli¬ 
cable  to  almost  all  marine  sediments. 

The  data  listed  in  the  tables  and  plotted  in  the  figures  do  not  include  attenuation 
in  calcareous  sediments.  There  is  a  small  amount  of  evidence  which  indicates  that  high- 
porosity  calcareous  "00/cs"  may  have  higher  attenuation  values  than  sediments  formed 
by  clay  minerals  of  the  same  general  porosity  and  grain  si/e.  This  might  In*  due  to  the 
nature  of  mterpattide  bonds  and  the  presence  oVtntaintni1eT.il  te>t>  in  the  calcareous 
oo/es. 

({fleets  of  Wetting 


When  a  dry  sand  is  wetted,  shear  velocity  and  rigidity  are  both  reduced  (Ref.  1‘). 

I  igs.  X  anil  d.  and  Ref.  2X).  and  the  log  decrement  increases  (Ref.  5b).  I’ilbeam  and 
Vaisnys  (Ref.  d4)  also  noted  that  lubrication  of  grains  in  a  granular  aggregate  generally 

increased  rvJ’diw  etterpy  TIv  roll?  til*  ripkUly  can  K*  see*  by  swMtttfVtafc • 

for  i;  in  I  q.  (4). 

=  A,(p)l;///j'2  (Id) 

Dynamic  rigidity  is  less  in  wet  sand  tlwu  in  dry  apparently  because  lubrication 
allows  grains  to  move  under  less  shear  stress  (Ref.  2K).  Rigidity  is  also  less  in  saturated 
porous  rock  than  in  dry  rock  (Ref. ‘>5).  p.  57X).  Hardin  and  Richart  noted  that  a  mere 
i  4  percent  iticnvMe  in  wrier  mulrul  can  reduce  rigidity  by  n  dnieM  m  1 5  pi  fee  IP  JA- 
perimenlally.  the  log  decrement  of  shear  waves  in  sand  and  rock  increases  upon  welting 


(Kefs.  5<>  and  95 ).  but  shear-wave  velocity  decreases  because  of  Ihe  decease  in  rigidity. 
Density  would  increase  with  wettint!.  Thus,  all  factors  in  liq.  ( I  3 1  result  in  increase  in 
a,  (alter  wettint!)  at  any  given  frequency.  A  variation  of  Iq.  (4).  \  =  indicates 

that  the  increase  in  the  log  decrement  upon  wetting  is  due  entirely  to  the  increase  in  as . 
because  I',  decreases.  The  apparent  reason  for  the  increase  in  attenuation  of  shear  waves 
upon  wetting  or  lubrication  is  that,  given  the  same  grain  si/e  and  porosity  (or  the  same 
number  of  interparticle  contacts),  the  reduced  rigidity  allows  greater  slippage  of  grains 
under  any  given  shear  stress,  resulting  in  greater  energy  losses  because  of  greater  inter¬ 
grain  friction. 

Effects  of  Pressure 

Although  attenuation  decreases  with  increasing  effective  pressure  (/.<*.,  overburden 
pressure),  as  discussed  below,  it  is  believed  that  the  laboratory  and  in  situ  measurements 
of  this  report  are  comparable  because  of  the  low  or  negligible  pressures  involved.  This 
was  also  a  conclusion  of  McCann  and  McCann  (Ref.  8)  after  comparing  their  in  situ  mea¬ 
surements,  at  depths  of  2  m  in  sediments,  with  the  laboratory  measurements  of  Shumway 
(Ref.  18).  The  in  situ  measurements  using  probes  (listed  in  Tables  I -4a)  were  at  sediment 
depths  from  about  0.3  to  2  m;  thus,  effective  pressures  were  very  small.  Laboratory 
samples  had  negligible  effective  pressures.  Consequently,  the  data  and  conclusions  of  this 
report  are  concerned  with  surlicial  sediments  under  very  low.  or  negligible  effective  pres¬ 
sures.  Some  of  the  measurements  at  seismic  frequencies  in  Table  4b  involved  greater 
effective  (overburden)  pressures,  but  these  were  not  involved  in  calculations  of  equations 
relating  attenuation  and  frequency. 

Increase  of  effective  pressure  increases  rigidity  and  shear-  and  compressional-wave 
velocity  in  rocks  (Refs.  59  and  9b)  and  sand  (Refs.  19  and  28);  but  in  both  rocks  and 
sand,  attenuation  and  log  decrement  decrease  (for  rocks,  see  Refs.  97.  98.  99;  for  sand. 
Refs.  27  and  5b).  As  various  authors  have  noted,  these  effects  are  apparently  due  to 
pressure  effects  on  grain  elastic  moduli  and  grains  in  harrier  contact,  a  situation  that  offers 
greater  resistance  to  shear  stress  (shear  modulus,  or  rigidity,  greater),  but  allows  less  inter¬ 
grain  movement,  which  constraint  reduces  energy  loss  through  intergrain  friction  (attenua¬ 
tion  less). 


PREDICTION  01  COMPRESSION  A  L-WAVE  ATTENUATION 

I  stimation  or  prediction  of  compressional-wave  attenuation  in  saturated  sediments 
is  uf  eHttshkwWe  importance  m  geophysics  ami  underwater  acoustics.  The  relationships 
between  attenuation  (as  expressed  by  the  constant  A)  and  grain  si/e  and  porosity  pre¬ 
viously  discussed  and  illustrated  d  igs.  3,  4.  and  5)  afford  a  very  simple  method  for  pre¬ 
diction  of  approximate  values  of  attenuation  at  most  frequencies  of  interest  in  geophysics 
and  underwater  acoustics  in  all  major  sediment  types  (with  the  possible  exception  of  cal¬ 
careous  sediments). 


In  this  report,  attenuation  is  related  to  frequency  in  the  form 


a  =  kf" 


(I) 


where 


a  is  attenuation.  dB  m 
I  is  frequency,  kll/ 

The  case  lias  been  made  that  the  exponent  of  frequency.  //.  is  close  to  I.  If//  is  taken 
as  I.  the  only  variable  in  the  equation  is  the  constant  A.  This  constant  varies  with  mean 
grain  si/e  d  igs.  3  and  4)  or  porosity  (Fig.  5). 

Mean  grain  si/e  of  a  sediment  can  be  easily  determined  in  both  wet  and  dry  sedi¬ 
ment;  porosity  can  be  determined  in  saturated  materials  In  the  absence  of  measurements 
both  properties  can  be  predicted,  given  general  environment  and  predicted  sediment  type. 
References  3  and  100  discussed  and  listed  averages  of  these  properties  for  various  sedi¬ 
ment  types  in  the  major  environments  and  discussed  methods  of  prediction.  Technically, 
porosity  should  be  a  better  index  of  attenuation  because  it  is  a  better  measure  of  the 
number  of  interparticle  contacts  (at  the  same  grain  si/e.  the  degree  of  packing  determines 
porosity  and  the  number  of  interparticle  contacts). 

I’he  method  for  estimating  or  predicting  approximate  values  of  attenuation  (given 
mean  grain  size  or  porosity )  follows;  a  numerical  example  is  included: 

( I )  Determine  or  predict  mean  grain  si/e  or  porosity  of  the  sediment, 
example:  mean  grain  si/o.  M:  ~  7.5  <p. 

(3)  Refer  to  the  mean  grain  si/e  or  porosity  versus  A  diagram  (Figs.  3.  4.  5).  and 
determine  the  value  of  A.  Regression  equations  for  these  data  are  listed  in 
the  figure  captions. 

Fxample:  Fig.  3  (mean  grain  si/e  versus  A  diagram)  yields:  A  =  0.07. 

(3)  Insert  A  into  Fq.  ( I ).  reproduced  above. 

Fxample:  insertion  of  the  value  of  A  yields:  n  =  0.07/ 1 

(4)  Compute  attenuation  (in  dlCm)  at  the  desired  frequency  (kll/) 

The  experimental  evidence  indicates  that  the  exponent  of  frequency  probably 
varies  in  saturated  sediments  between  0.0  and  1.1.  Theoretically,  or  statistically,  the 
value  might  be  slightly  less  than  I  (Refs.  45  and  51).  Use  of  a  value  of  A.  assuming  n 
is  I  (rather  than  0.0  or  I.I).  results  in  small  differences  in  computed  attenuation  at  lower 
frequencies,  those  of  most  interest  in  geophysics  and  underwater  acoustics.  For  example. 


al  .3  kll/.  if  0.07  is  taken  as  k.  ami  n  is  varied  between  ().*>.  1.0.  and  l  .l.  the  computed 
attenuations  are  0.l*>.  0.21.  and  0.2.3  dli  m:  at  I  -  50  II/.  the  values  of  attenuation  are 
0.005.  0.004.  and  0.00.3  dB  m.  In  a  fine  sand,  if  k  is  0.5.  the  values  of  attenuation  at 
.3  kll/  (and  n  =  ().‘>.  1.0.  and  l.l )  are  1  . 34,  1.50.  and  1  .(>7  dB  m:  at  50  II/.  attenuations 
are  0.0.34.  0.025.  and  0.01*)  dB  m. 
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